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ABSTRACT 


Due  to  inherent  turbulence,  the  atmosphere  has  a  temporally  and 
spatially  variable  refractive  index,  which  degrades  propagating 

2  'V 

electromagnetic  radiation.  is  a  key  parameter  for  describing 

refractive  variations.  Data  was  analyzed  from  several  instruments 

involved  in  an  electro-optics/meteorology  experiment:  at  Penn  State 

University!  -,  scintillometer  (RADC,  AFGL,  NPS)  ,.  thermosonde  (AFGL)  , 

radar  vertical  profiler  (ESU) ,  and  instrumented  aircraft  (ARA) . 

2  2 

The  aircraft  measured  C^,  and  using  hot-  and  cold-wire 

"  / 

sensors  and  FM  recording  apparatus.  The  taped  data  was  processed  via 

FFT  to  produce  one -dimensional  variance  spectra  (wavenumber  range  0.01 

/ 

to  10  m  ) .  Flights  usually  produced  a  10  km  vertical  profile 

processed  to  give  roughly  0.5  km  resolution.  Spectral  editing  was 

based  on  regression  analysis  fit  to  -5/3  frequency  dependence.  A 

2 

majority  of  spectra  showed  the  classic  inertial  subrange.  Cn  and 

2  2 

e  were  calculated  from  and  Cu  ,  respectively.  In  active 
regions,  the  following  relationship  can  be  derived  (Fairall  and 
Markson  1984) : 

C T  /  =  1.6  RL/  (Pr  -  Rt)  (*/g)  (ae/dz) 

2  2 

A  scatter  plot  of  C^,  vs.  ( 6/g  86 /d z)  in  the  free 

troposphere  showed  a  range  of  values  from  0.3  to  10  for  high  and  low 

2  2  2 
values  of  C^,  ,  respectively.  versus  showed  high 

2  2 

correlation,  with  log  “1.7  log  C,p  .  Plots  of 

/Q^ (6 /g  ae/dz)  vs.  e/(v  N)  (Gregg  1987)  approached  an 
ordinate  value  of  0.4  at  high  activity  levels.  Values  of  N  and  S  were 


2 

used  to  obtain  values  of  C  and  t  with  the  model  of  VanZandt  et 

n 

2 

al.  (1981).  Comparisons  of  Cn  profiles  showed  average  ratios  of 

each  to  the  thermosonde  profile  were:  aircraft  8,  model  6.8,  profile 

8.1,  and  the  scintillometer  1.1.  Numerical  estimates  of  8  from 

o 

2 

Cn  profiles  differed  from  the  aircraft  by  an  average  of  17%  for 

the  thermosonde  and  scintillometer.  Estimates  of  r  differed  by  35% 

o  J 

for  the  thermosonde.  Aircraft-measured  and  model-predicted  values  of 
e  differed  by  several  orders  of  magnitude. 
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1  ABSTRACT 


|  Due  to  its  inherently  turbulent  nature,  the  atmosphere  has  a 

|  temporally  and  spatially  variable  refractive  index,  which  degrades 

! 

|  propagating  electromagnetic  radiation.  The  refractive  index  turbulent 

I  2 

f  structure  constant,  C  ,  is  a  key  parameter  for  describing 

I  n 

i  « 

,  refractive  variations.  C  ‘  can  be  inferred  from  optical 

j  turbulence  (scintillometer) ,  temperature  turbulence  (aircraft 

|  instrumentation  or  thermosonde)  or  radar  backscatter  (profiler) 

!  2 

,  measuren ents .  Vertical  integrals  of  C  give  transverse  coherence 

i 

!  length,  rQ,  and  isoplanatic  angle,  8q.  These  can  also  be  measured 


remotely  by  an  rQ  scintillometer  and  isoplanometer ,  respectively. 

During  an  atmospheric  optics/meteorology  experiment  (acronym 

EWAK)  conducted  at  Penn  State  University  primarily  during  April  and 

May  of  1986,  data  was  collected  by  all  of  the  aforementioned 

instruments.  An  instrumented  research  aircraft  was  used  to  measure 

2 

vertical  profiles  of  temperature  and  velocity  turbulence  (C^  and 
2 

)  and  other  meteorological  variables.  The  turbulence 
instrumentation  aboard  the  aircraft  consisted  of  cold-wire  and 
hot-wire  sensors  and  FM  recording  apparatus.  The  taped  data  was 
processed  via  FFT  to  produce  one-dimensional  variance  spectra 
(wavenumber  range  0.01  to  10  m  b .  Flights  usually  produced  a  10  km 
vertical  profile;  the  data  was  processed  to  give  roughly  0.5  km 
vertical  resolution  (similar  to  that  of  the  radar  wind  profiler) . 

Editing  was  based  on  percent  error  between  a  regression  analysis 
and  theoretical  -5/3  frequency  dependence  of  the  spectra.  A  majority 
of  spectra  showed  good  evidence  of  the  classic  inertial  subrange  with 


-5/3  slope.  C^  and  were  calculated  from  the  regression 

fit  to  the  inertial  subrange  power  spectral  density.  The  rate  of 

dissipation  of  turbulent  kinetic  energy,  e,  was  calculated  using  the 

2 

so-called  Corrsin  relation  for  velocity  turbulence.  Optical 

2 

was  calculated  from  . 

Considerable  interest  has  developed  in  models  that  relate 

micro turbulence  parameters  to  the  mean  gradients.  One  such  model, 

proposed  by  VanZandt  et  al .  (1978;  1981),  has  as  a  key  variable  the 

gradient  Richardson  number,  .  In  actively  turbulent  regions, 

2  2 

can  be  related  to  C^,  and  as  follows: 

c\  /  1.6  Rj/  (Pr  -  R.)  (tf/g)  {de/d z) 

2  ? 

A  scatter  plot  of  CT  versus  (9/ g  d9/d z)  Cu  in  the  free 

troposphere  showed  a  range  of  values  from  0.3  to  10  for  high  and  low 

2  2 

values  of  CT  ,  respectively.  Scatter  plots  of  CT  versus 

2 

Cu  showed  high  correlation  of  these  two  parameters,  with  log 
2  2 

-  1.7  log  CT  .  R^-N/S  and  plots  of  N  versus  S  roughly 

corresponded  to  the  curve  depicted  by  Fairall  and  Markson  (1985)  . 

2  2 

Plots  of  C^  /  (0/g  dd/dz)  versus  a  turbulent  activity 

parameter  (Gregg  1987)  clearly  showed  agreement  with  the  suggested 
activity  levels  and  associated  values.  At  high  activity  levels,  the 
ordinate  value  approached  0.4  (equivalent  to  1.6  when 
approximating  P  as  1  and  R^  as  0.25). 

2 

Values  of  N  and  S  were  used  to  obtain  values  of  C  and  « 

n 

2 

with  the  model  of  VanZandt  et  al .  Comparisons  of  C^  profiles 
measured  by  the  various  instruments  showed  good  agreement  between 
scintillometer  and  thermosonde,  and  between  the  VanZandt  et  al .  model 


and  the  aircraft  data.  The  profiler  (converted  optical  values)  had 
some  agreement  with  the  thermosonde.  Average  ratios  of  each  profile 
to  the  thermosonde  profile  were:  aircraft  8,  model  6.8,  profiler  8.1, 
and  the  scintillometer  1.1. 

2 

Numerical  integration  of  gave  estimates  of  rQ  and  8q. 

Values  from  scintillometer,  thermosonde  and  aircraft  did  not  show 

clear  agreement.  Based  on  the  aircraft  values,  the  scintillometer  and 

thermosonde  6  values  differed  by  an  average  of  17%.  The 

thermosonde  rQ  values  differed  by  an  average  of  35% . 

2 

Despite  good  agreement  on  ,  the  aircraft  data  and  the 
VanZandt  et  al.  model  output  clearly  disagreed  on  vertical  profiles  of 
£,  in  some  cases  by  several  orders  of  magnitude. 


Overall,  EWAK  provided  a  chance  to  compare  methods  of  measuring 

2 

Cn  ,  as  well  as  an  opportunity  to  obtain  microturbulence  data  for 
examining  relations  to  the  mean  gradient  structure. 


oDI 

I 

& 


% 


I 

1 

'V 


$ 


k’A 


TABLE  OF  CONTENTS 


ABSTRACT . 


LIST  OF  TABLES . viii 


LIST  OF  FIGURES. 


ACKNOWLEDGEMENTS . 


Chapter 


INTRODUCTION . 


BACKGROUND  ON  ATMOSPHERIC  PROPAGATION. 


Refractive  Index. 


2.2  Propagation  Effects. 


2.3  Propagation  Parameters  and  C 


Measurements . 


MICROTURBULENCE  THEORY. 


The  Refractive  Index  Structure  Function. . .  10 


C  From  Radar  Measurements . 
n 


2  2 

C^,  and  Cu  Relationships .  13 

The  Turbulent  Kinetic  Energy  Dissipation 
Rate .  16 


Turbulence  Size  Scales. 


THE  EWAK  EXPERIMENT. 


Overview . 


Aircraft  Instrumentation. 


Additional  Instrumentation. 


5 


DATA  PROCESSING 


25 


TABLE  OF  CONTENTS 
(continued) 


Chapter 

6  ANALYSIS  METHODS . 

6.1  Aircraft  Turbulence  Data  Analysis.... 

6.2  Mean  Aircraft  Meteorological  Profiles 

6.3  Other  Data  Sources . 


7  RESULTS  AND  DISCUSSION . 

2  2 

7.1  Analysis  of  C^,  /C^  . 

7.2  Mean  Gradients  and  Richardson  Number. 

7 . 3  Profiles  of  C  ^ . 

n 

7.4  Propagation  Parameters . 

7.5  VanZandt  et  al .  Model  for  € . 

8  SUMMARY  AND  CONCLUSIONS . 

REFERENCES . 


■■  -.  V  (.■  •v- 


LIST  OF  TABLES 


Table 


Aircraft  flight  logistics  and  corresponding  data 
from  other  instruments  during  EWAK . 


Representative  editing  results  for  each  flight  and 
in  total . 

Estimates  of  the  height  of  the  PBL  based  on  the 
strength  of  the  signal  - to-noise  ratio  exhibited  in 
the  spectra . 


Calculated  vertical  profiles  of  Richardson  number 
based  on  R.-N/S  and  using  N  and  S  values  shown  in 
Figure  17  .  ^ .  120 

The  arithmetically  averaged  Richardson  number  at 
each  height  for  all  vertical  flights .  123 


X  vertical  profiles  for  each  flight  during  which 
the  profiler  operated,  calculated  at  each  aircraft 
observation  level  in  MSL  altitude  with  aircraft- 
measured  temperature  and  specific  humidity 
profiles .  129 


X  vertical  profiles  for  each  flight  during  which 
the  profiler  operated,  calculated  at  each  profiler 
observation  level  in  MSL  altitude  with  aircraft- 
measured  temperature  and  specific  humidity 
profiles .  132 


Arithmetically  averaged  X^  vertical  profiles  from 
the  Table  7  profiles .  135 

AFGL-calculated  optical  parameters  from  thermosonde 
data .  147 


Average  of  AFGL-calculated  optical  parameters  from 
scintillometer  data .  148 


11  Optical  parameters  from  aircraft  data  calculated 
with  the  Murphy  and  Battles  assumptions . 


149 


LIST  OF  FIGURES 


Ml 


Figure 


Voltage  variance  spectra  from  first  FFT  program 

(128  points) .  31 

Voltage  variance  spectra  produced  by  the  signal 
analyzer's  1600  point  FFT .  34 


The  measured  2  Hz  high  pass  temperature  turbulence 
channel  filter . 

The  measured  5  Hz  high  pass  velocity  turbulence 
channel  filter . 

The  derived  atmospheric  2  Hz  high  pass  temperature 
turbulence  channel  filter . 


The  derived  atmospheric  5  Hz  high  pass  velocity 
turbulence  channel  filter . 


Voltage  variance  spectra  showing  editing  procedure.  41 

Voltage  variance  spectra  showing  editing  procedure.  44 

Voltage  variance  spectra  showing  editing  procedure.  47 

Potential  temperature  (theta)  vertical  profiles  as 
measured  by  aircraft  instrumentation .  55 

Specific  humidity  vertical  profiles  as  measured  by 
aircraft  instrumentation .  60 


Mean  wind  speed  (solid  line)  and  direction  (broken 
line)  as  measured  by  the  profiler . 


AFGL  thermosonde-measured  potential  temperature 
vertical  profiles  correlating  to  aircraft  flights..  66 

A.FGL  thermo  sonde -measured  specific  humidity 

vertical  profiles  correlating  to  aircraft  flights..  67 

A.FGL  thermosonde -measured  vertical  profiles  of  mean 

wind  speed  (solid  line)  and  direction  (broken 

line) .  68 

AFGL  t^ermosonde - calculated  vertical^prof iles  of 

raw  -optical  based  on  measured  C^,  .  69 


Figure 


LIST  OF  FIGURES 
(continued) 


AF^L  scintillometer-measured  vertical  profiles  of 

C  -optical  taken  during  the  aircraft  flight 

time  period .  71 

RA^C  scintillometer-measured  vertical  profiles  of 

C  -optical  taken  during  the  aircraft  flight 

time  period .  72 

Vertical  profiles  of  C  ^  ,  C_^  ,  and  C  ^ .  74 

Y  n  T  ’  u 


A  comparison  for  each  vertical  flight  of  C,^ 


versus  Cu  for  spectra  above  the  PBL. 


A  comparison  for  each  level  flight  of 
versus  for  spectra  above  the  PBL. . . . 


2  2 

A  composite  comparison  of  C^,  versus 
for  all  vertical  profiles . 


2  2 

A  composite  comparison  of  C^  versus 
for  all  level  profiles  . 


2  2 

A  composite  comparison  of  C^  versus  C 

normalized  by  the  expression  [T/g  (5T/$z  +  T) ] .  104 

2 

Comparison  fo£  each  vertical  profile  of  C^, 

divided  by  C  normalized  as  in  Figure  20 

versus  a  turbulent  activity  parameter,  e/(u  N) .  106 


2  2 

A  composite  comparison  of  C^,  divided  by  Cu 
normalized  as  in  Figure  20  versus  a  turbulent 
activity  parameter  s/(i/  N)  for  all  vertical 
profiles .  115 

Comparisons  of  the  VanZandt  et  al .  parameters  N, 

the  Brunt-Vaisala  frequency  squared,  and  S,  the 

shear  squared .  116 

A  composite  comparison  of  N  versus  S  for  all 

flights  shown  in  Figure  23 .  119 

Comparison  of  the  VanZandt  et  al.  model  vertical 

profile  prediction  of  C  -optical  (solid 

line)  and  the^aircraft  -derived  vertical 

profile  of  C^  -optical  (broken  line) .  124 


124 


M*mU\ 


.1^1 


Figure 


LIST  OF  FIGURES 
(continued) 


A  comparison  of  C  -radar  as  derived  from 
profiler  (solid  line)  and  from  aircraft 
-optical,  multiplied  by  the  conversion 
factor  X  to  give  C  -radar  (broken 
line) . n . 


Comparison  of  vertical  profiles  of  -optical .  141 

2 

Comparison  of  vertical  profiles  of  C^  -optical .  143 

Comparison  of  vertical  profiles  of  epsilon  (e),  the 
turbulent  kinetic  energy  dissipation  rate, 
predicted  by  the  VanZandt  et  al .  model  (solid 
line)  and  measured  by  the  aircraft  (broken  line) . . .  151 


Average  vertical  profile  of  eta  ( rj )  ,  the 
turbulence  inner  size  scale,  from  aircraft- 
measured  epsilon  profiles . 


•*.W 


V/AAV./J 


V.V.VV 


ACKNOWLEDGEMENTS 


First  and  foremost,  I  would  like  to  express  my  sincerest 
gratitude  to  my  thesis  advisor,  Dr.  Christopher  W.  Fairall.  Without 
his  gentle  guidance  and  patient  corrections,  this  work  could  not  have 
been  completed.  I  also  acknowledge  Dr.  Dennis  W.  Thomson  for  his 
helpful  suggestions  concerning  my  work.  A  special  thanks  to  my  fellow 
graduate  students,  Mark  Miller  and  James  Edson,  and  their  computer 
expertise  that  they  so  willingly  shared. 

With  regard  to  my  program  coursework,  I  thank  Dr.  Craig  F.  Bohren 
and  his  humor  for  putting  a  smile  on  my  face  whenever  I  needed  one  and 
for  always  reminding  me  of  the  shear  joy  of  studying  a  physical 
science . 

Thanks  also  to  all  those  involved  in  EWAK  who  provided  their 
time,  expertise,  instruments,  and  data,  and  who  answered  my 
multitudinous  questions:  Dr.  Ralph  Markson  (ARA) ,  Capt  Mike  Moss 
(Florida  State  University),  Mr.  Donald  Stebbins  (RADC) ,  Mr.  Paul  Cook 
(RADC) ,  Mr.  James  Brown  (AFGL) ,  Mr.  Robert  Beland  (AFGL) ,  Mr.  Edmund 
Murphy  (AFGL),  Dr.  Donald  Walters  (NPS)  and  Capt  Jack  Davidson  (AFWL) . 

Finally,  I  would  like  to  thank  my  loving  husband,  Jeff  Beecher, 
for  providing  "R  and  R"  when  needed,  as  well  as  always  willingly 
acting  as  my  sounding  board  for  ideas  and  my  source  of  practical 
mechanical  knowledge. 

This  research  was  sponsored  in  part  by  the  U.  S.  Air  Force  under 
contract  #AFOSR- 86 -0049 . 


CHAPTER  1 


INTRODUCTION 

For  many  reasons,  there  is  considerable  interest  in  the  vertical 

2 

profile  of  the  refractive  index  structure  parameter,  C^  .  This 
parameter  quantifies  effects  of  refractive  index  fluctuations 
important  for  electromagnetic  radiation  systems  whose  signals 
propagate  through  the  turbulent  atmosphere.  Examples  of  such  systems 
are  satellite  communications,  remote  sensing  such  as  clear  air  radars, 
astronomical  observation,  and  certain  weapons.  In  situ  measurements 
are  a  necessary  first  step  to  set  up  a  database  for  verifying  physical 
and  numerical  models.  Both  are  required  to  facilitate  the 
interpretation  of  remote  sensing  measurements  and  to  achieve  the  goal 
of  prediction  based  on  readily  available  atmospheric  data  and  models . 

This  project  was  an  analysis  of  the  first  segment  of  the  EWAK 
(E  xperiment  W  ithout  AK  ronym)  electro-optical/meteorological 
experiment  at  Penn  State  University.  The  purpose  of  EWAK,  conducted 
principally  in  April  and  May  of  1986,  was  to  simultaneously  operate  at 
the  same  location  several  different  instruments  to  measure  the 
refractive  index  structure  parameter,  both  in  situ  and  remotely, 
throughout  different  synoptic  and  diurnal  conditions. 

Specifically,  this  thesis  describes  a  detailed  analysis  of 
research  aircraft  measurements  of  temperature  and  velocity  turbulence 
taken  with  the  instrumented  Beechcraft  Baron  of  the  Airborne  Research 
Associates  of  Weston,  Massachusetts.  From  these  two  measurements  can 

be  derived  the  temperature  and  velocity  structure  parameters, 

2  2  2 
and  C  .  Optical  C  is  a  function  of  C„  . 
u  r  n  1 


The  turbulent 


kinetic  energy  dissipation  rate,  e,  is  related  to  .  Another 

2  2 

important  relationship  examined  was  that  between  and  Cu  . 


CHAPTER  2 


BACKGROUND  ON  ATMOSPHERIC  PROPAGATION 

2.1  Refractive  Index 

A  turbulent  fluid  creates  three  general  effects:  it  imposes 
varying  forces  on  any  objects  embedded  in  the  fluid  or  in  the  fluid 
path,  it  generates  fluxes  of  fluid  properties  (e.g.,  temperature  or 
momentum),  and  it  creates  (by  density  fluxes)  variations  in  the 
refractive  index  of  the  fluid  (Panofsky  and  Dutton  1984) .  These 
refractive  index  variations  will  affect  electromagnetic  (EM)  and  sonic 
radiation  propagating  through  the  fluid.  The  atmosphere,  being  a 
turbulent  fluid,  has  a  temporally  and  spatially  variable  refractive 
index  structure.  The  atmosphere's  turbulent  refractive  index 
fluctuations  have  plagued  astronomers  since  some  of  the  earliest 
recorded  observations.  The  twinkling  of  stars  is  one  of  the  most 
obvious  visible  manifestations.  Today,  our  atmospheric  uses  of  EM 
radiation  are  far  more  diverse  than  merely  receiving  visible  starlight 
with  our  eyes  or  at  a  telescopic  aperture.  Hence,  the  need  exists  to 
characterize  and  predict  the  refractive  index  and  its  variations  over 
arbitrary  pathlengths  and  view  angles. 

Henceforth,  in  this  thesis,  unless  it  is  specifically  otherwise 
stated,  the  discussion  will  be  concerned  with  propagating 
electromagnetic  radiation.  Refraction  can  be  considered  as  a  type  of 
scattering  (i.e.,  deviation  from  the  forward  direction).  Turbulent, 
non-homogeneous  variations  of  the  atmospheric  refractive  index  are  on 
the  order  of  one  part  in  10^.  This  results  in  very  small  scattering 
angles,  for  which  depolarization  and  intensity  attenuation  generally 


can  be  neglected  (Strohbehn  1978).  The  real  part  of  the  complex 
refractive  index,  n,  of  a  medium  can  be  thought  of  as  a  measure  of  the 
medium's  effect  on  wave  propagation  as  compared  to  propagation  in  a 
vacuum:  n(medium)=c (vacuum) /c (medium) ,  where  c  is  the  phase  velocity 

of  the  wave.  Refractivity ,  usually  denoted  by  N,  is  (n-1)  X  10^. 

The  phase  speed  of  the  wave  is  altered  (slowed  in  air  compared  to 
vacuum,  thus  n>l) ,  but  the  energy  E  and  frequency  u  remain  the  same, 
E~hi/  (Hecht  1986) .  Using  geometric  optics  theory,  this  can  be  shown  to 
result  in  a  deviation  of  the  beam's  path  which  is  referred  to  as  the 
phenomenon  of  refraction  (in  molecular  dipole  theory,  it  is  referred 
to  as  a  scattering  event) . 

Two  factors  determine  the  phase  speed  (and  thus  refractive  index) 
in  a  medium:  density  (number  of  atoms,  molecules  or  particles  per  unit 
volume)  and,  since  each  molecule  has  a  wavelength  dependent 
polarizability,  molecular  composition  (Hecht  1986).  At  optical  or 
shorter  wavelengths,  temperature  variations  effectively  dominate 
refractivity  changes  since  atmospheric  molecules  have  low  rotational 
polarizability  at  these  wavelengths  (Balsley  and  Gage  1980).  However, 
at  longer  wavelengths,  specifically  microwave  radar  wavelengths,  both 
temperature  and  the  relative  concentration  of  various  molecular 
species  with  permanent  dipole  moments  determine  the  refractive  index 
gradients  affecting  radiation  transmittance  (Wesely  1976;  Balsley  and 
Gage  1980).  In  the  boundary  layer,  water  vapor  content  is  often 
highly  variable,  and  thus,  may  dominate  refractivity  changes  (Balsley 
and  Gage  1980) . 

If  temperature  and  molecular  composition  varied  negligibly  from 
the  standard  atmosphere,  then  the  refractive  index  profile  could 


easily  be  calculated.  The  importance  of  turbulence  is  its  resultant 
mixing  (Strohbehn  1978).  Mixing  due  to  turbulence  causes  the 
variations  in  space  and  time  of  the  refractive  index  profile.  It  is 
the  random  nature  of  turbulence  that  introduces  chaos  into  the  density 
profile . 

2.2  Propagation  Effects 

Variations  in  refractive  index  can  cause  at  least  six  types  of  EM 
propagation  degradation:  beam  steering,  image  dancing,  beam  spread, 
spatial  coherence  degradation,  temporal  coherence  degradation,  and 
scintillation  (Dewan  1980) . 

Beam  steering  is  deviation  from  the  line-of-sight  (i.e.,  the  most 
direct  path  from  source  to  receiver) .  The  transverse  deviation  of  a 
narrow  beam  increases  with  subsequent  distance  from  the  deviating 
medium.  If  a  turbulent  eddy  is  larger  than  the  width  of Vhe  beam,  it 
may  act  as  a  "prism."  This  effect  can  be  quite  significant  for 
earth -to -space  uplinks,  because  the  planetary  boundary  layer  (PBL)  is 
highly  turbulent  (Dewan  1980) . 

Image  dancing  is  the  modulation  of  the  arrival  angle  of  the  whole 
beam  wavefront.  It  can  cause  the  focal  plane  for  the  beam  to  be 
angled  to  a  receiver  normal  to  the  beam,  causing  photographic  blurring 
or  necessitating  larger  apertures  with  lower  signal- to-noise  ratios 
(Dewan  1980) . 

Beam  spreading  is  caused  by  scattering  at  small  angles  near  the 
edge  of  the  beam.  This  occurs  when  the  beamwidth  is  greater  than  the 
turbulent  eddies.  This  phenomenon  occurs  most  on  the  space  -  to  -  earth 


downlink  (Dewan  1980) . 


■BBS iff! 


At  large  distances,  a  spherical  wave  from  a  point  source  can  be 
approximated  as  a  plane  wave.  According  to  Huygen' s  model,  a  single 
wavefront  can  be  thought  to  consist  of  an  infinite  number  of  points 
each  emitting  a  "wavelet."  The  wavelets  add  together  to  give  the 
total  wave.  If  the  wavelets  have  a  constant  phase  relation,  the 
wavefront  is  said  to  be  coherent.  When  phase  difference  randomly  and 
continuously  changes,  the  wavefront  becomes  incoherent  (Hecht  1986). 
When  a  pathlength  difference  is  introduced  across  a  coherent  wavefront 
(e.g.,  refraction),  a  Fraunhoffer  diffraction  pattern  may  be  seen  in 
which  maxima/minima  occur  where  the  amplitudes  of  wavelets 
contructively/destructively  interfere.  Thus  diffraction  is  one  method 
of  testing  coherence  (Hecht  1986) . 

A  coherent  source  can  be  viewed  as  emitting  "pulses,"  or  wave 
packets  (this  model  combines  wave  and  particle  concepts).  No  source 
is  perfectly  monochromatic,  thus  each  wave  pulse  has  a  frequency 
bandwidth.  The  inverse  of  the  pulse  frequency  band  is  the  "period" 
band  or  coherence  period  (Hecht  1986).  Temporal  (or  longitudinal) 
coherence  occurs  within  the  coherence  period.  Spatial  (or  horizontal) 
coherence  occurs  across  a  wave  pulse  emitted  from  a  coherent  source. 
Spatial  degradation  is  loss  of  phase  coherence  across  a  wavefront, 
causing  spot  blurring  on  photographic  images.  It  also  distorts 
photomixing  processes  which  rely  on  mixing  of  an  incoming  to  a  local 
signal  to  form  a  "beat."  Temporal  degradation  is  a  loss  of  phase 
coherence  during  the  coherence  period,  thus  distorting  amplitude 
modulated  signals  (Dewan  1980) . 

Scintillation  is  a  coherence  effect  describing  amplitude  changes 


at  a  focal  plane  (interference  patterns).  Spatial  scintillation  being 


fluctuations  across  a  beamwidth  (or  wavefront) ,  temporal  scintillation 
being  fluctuations  in  time  of  the  total  beam  amplitude  (signal  fading) 
(Dewan  1980). 

2 

2.3  Propagation  Parameters  and  Measurements 


A  variety  of  theories  have  been  developed  for  these  degradation 

effects.  In  general,  an  atmospheric  turbulence  parameter,  the 

2 

refractive  index  structure  constant,  ,  must  be  specified  along 

the  propagation  path.  By  measuring  some  relevant  propagation 

parameter  (e.g.,  scintillation)  the  theory  can  be  "inverted"  to 
2 

determine  C 

n 

One  such  instrument,  the  NOAA  (National  Oceanic  and  Atmospheric 

2 

Administration)  scintillometer,  measures  Cn  -optical  remotely, 

based  on  starlight  amplitude  fluctuations.  As  one  of  the  six 

distortions  described  previously,  scintillation  has  been  quantified 

2 

using  an  integral  of  Cn  (z) ,  where  z  is  the  altitude.  Tatarskii 
has  shown  for  temporal  scintillations,  time  averaged  amplitude 
fluctuations  (Dewan  1980) ,  that 


< [ ln(A/A  )]2>  -  0.56  k7/6 


„2 ,  .  5/6, 

C  ( z )  z  '  dz 
nv  ' 


(2.1) 


where  A  is  the  signal  amplitude,  Aq  is  the  mean  amplitude,  k  is  the 

optical  wave  number,  and  the  angle  brackets  denote  a  time  average.  A 

2 

stellar  scintillometer  measures  (d.),  the  normalized  variance 

of  the  spatially  filtered  signal  irradiance  for  3  spatial  wavelengths, 

2 

d^  .  Then  a  profile  of  is  calculated  for  seven  predetermined 

levels,  based  on  seven  weighting  functions  (only  four  of  which  are 


independent) 


(2.2) 


C^-BE  R.  K.8/3  a\  (cL) 
n  ,,ii  11 
l-l 


where  R^  weight  is  given  to  the  ith  spatial  wavenumber 
(K^-2?r/cL)  and 


B  -  1.87  X  10'13  [  Wc(z)  dz  l’1 

o 


(2.3) 


where  Wc(z)  is  a  composite  path  weighting  function  (Ochs  et  al . 


1976).  The  seven  predetermined  levels  are  line-of-sight  distance  from 


the  instrument,  so  that  actual  height  above  ground  depends  on  the 


zenith  angle  of  the  observed  star.  Some  improvements  to  the  NOAA 


scintillometer  were  made  by  AFGL  (Air  Force  Geophysical  Lab)  to  their 


instrument.  They  are  described  in  more  detail  by  Murphy  and  Battles 


(1986)  . 


An  NPS  (Naval  Postgraduate  School)  rQ  scintillometer  optically 


measures  the  transverse  coherence  length,  rQ,  of  the  image  in  the 


focal  plane  (Donald  Walters,  telephone  conversation,  October  1987). 


rQ  is  a  scalar  measure  (in  units  of  aperture  size)  of  spatial 


wavefront  coherence.  Crudely  stated,  the  resolution  of  an  image 


cannot  be  improved  by  increasing  a  single  aperture  beyond  rQ  (Dewan 


1980) .  This  helps  delineate  the  aperture  size  for  adaptive  optics 


elements.  It  has  been  shown  that 


rQ  =  18.5  A1'2  (  C2(z)  dz  )*3/5 

o 


(2.4) 


2  -2/3 

with  r  in  cm,  wavelength  (A)  in  m,  C  in  m  and  altitude 

o  °  n 


in  m.  Using  this  integral,  rQ  can  also  be  calculated  from  profiles 


of  C  (Fried  and  Mevers  1974) . 
n 


I 

ft 
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An  AFWL  (Air  Force  Weapons  Lab)  isoplanometer  optically  measures 


8  ,  the  isoplanatic  angle.  8  is  a  measure  of  angular  coherence 

for  an  object  with  angular  size  (vs.  a  point  source)  that  delineates 

the  region  through  which  the  object  can  be  viewed  coherently  without 

adaptive  optics  (Loos  and  Hogge  1979;  Murphy  and  Battles  1986).  8 

2 

can  also  be  calculated  from  a  vertical  profile  of  (Murphy  and 

Battles  )986) 


9  -  58000  A 

o 


1.2 


z5/3  C2 (z)  dz  f3/5 
n  1 


(2.5) 


with  8q  in  /xrad,  and  A,  z  and  in  the  same  units  as  for 


Notice  that  the  three  parameters  discussed  here  have  different 

2  5 /6 

'ropagation  path  weighting  functions  (a^  as  z  '  ,  rQ  as  1.0, 

5 /3 

9q  as  z  ).  Thus,  for  ground- to  -  space  optical  paths,  rQ  tends 

to  be  dominated  by  the  boundary  layer,  9q  by  the  stratosphere,  and 

2 

ctt  is  somewhere  in  between. 


CHAPTER  3 


MICROTURBULENCE  THEORY 

3.1  The  Refractive  Index  Structure  Function 

The  atmosphere  is  predominantly  hydrostatically  stable  above 
the  planetary  boundary  layer  (PBL) .  Nevertheless,  observations,  such 
as  contrails  and  rocket  trails,  have  clearly  shown  that  discreet 
turbulent  layers  exist  (VanZandt  et  al .  1981).  Although  turbulence 
and  byproducts  of  turbulence  (fossil  eddies)  are  undoubtedly  present, 
they  exist  only  on  an  intermittent  basis. 

while  atmospheric  turbulence  is  inherently  anisotropic,  at  small 
size  scales,  where  viscous  dissipation  occurs,  turbulence  can  be 
considered  isotropic.  For  isotropic  and  homogeneous  turbulence,  the 
structure  function,  D,  depends  on  spatial  scalar  distance,  £,  and  is 
defined  as  the  mean  squared  difference  of  the  measured  passive 
parameter  (e.g.,  temperature)  at  two  points  in  space  (Hinze  1975). 

Dt<£>  -  <(T(r)  -  T(r+0)2>  (3.1) 

where  the  angle  brackets  indicate  an  averaged  quantity.  In 
anisotropic  turbulence,  D  also  depends  on  the  orientation  of  the 
separation  (Hinze  1975)  . 

The  autocorrelation  R  is  defined  as  the  average  of  the  product  of 
deviations  from  the  mean  (e.g.,  T')  at  two  points  in  space. 

Rt(?)  -  <T '  (r)  T'  (r+0>  (3.3) 

In  the  limit,  as  £  approaches  infinity,  the  structure  function 
approaches  twice  the  variance.  Variance  is  equal  to  the 


autocorrelation  function  at  £-0  (Tatarskii  1971) .  Thus 


2  dT(0  -  Rt(0)  -  R*r(0  (3.3) 

In  the  inertial  subrange,  it  can  be  shown  that  the  structure 
function  has  a  power  law  relation  to  separation  with  a  constant 
designated  (e.g.,  C^2)  (Hinze  1975). 

Dt<£>  -  C2  |2/3  (3.4) 

Similarly,  we  can  define  for  refractive  index, 

0n(O  -  C2  £2/3  (3.5) 

The  structure  constant  is  only  a  "constant"  in  that  it  is  independen 

of  the  separation  distance;  it  is  still  variable  over  time  and  space 

2 

[e.g.,  CT  (x,y,z,t)].  The  optical  refractive  index  structure 
constant  can  be  related  to  the  temperature  and  humidity  structure 
constants  (Nastrom  et  al .  1981). 

*  hh  ^ 

where  f^  =  0.60  P/T2,  and 

f2=  -  1.2  Pq/T2  6  -  77.6E-6  P/  T  8  (3.7) 

with  pressure  in  mb , and  where  q  is  specific  humidity.  Neglecting  th 
effects  of  humidity,  this  reduces  to  (Tatarskii  1971) 

C2  (optical)  -  (79  X  10~6  (P/T2))2  C2  (3.8) 

2 

During  the  EWAK  experiment,  the  aircraft -measured  allowed 

2 

computation  of  "dry"  C  -optical.  A  Lyman-alpha  detector  on  the 


with  (Fairall  and  Thomson  1985) 


X2  -  [1-  7750  q'  (q/T)  (g/N)  +  15500  q/T]2  (3.13) 

2  2 

3.3  C^,  and  Relationships 

Dimensional  analysis  shows  that  in  the  inertial  subrange  of 

-  5  /  3 

isotropic  turbulence,  the  variance  spectral  density  obeys  a  k  / 
wavenumber  dependence . 

^T(kx)  -  /3t  xt  *'1/3  k‘5/3  (3.14) 


where  is  an  empirical  constant,  is  the  rate  of  dissipation 
of  one  half  the  temperature  variance,  e  is  the  dissipation  rate  of  th 
turbulent  kinetic  energy,  and  k^  is  the  spatial  wavenumber  (k^  - 
2n  X  ■*■)  (Panofsky  and  Dutton  1984)  .  Since  the  one  dimensional 
variance  spectrum  is  defined  as  the  Fourier  transform  of  the 
autocorrelation  (Panofsky  and  Dutton  1984) , 


*  T 


CO 

RT(f1)  cos  (k1f1)  df 


(3.15) 


Jo 

it  can  be  shown  using  the  relation  of  the  structure  function  to 
variance  and  the  autocorrelation,  that 


<t>T( kx)  =  0.25  C2  k^5/3 


(3.16) 


in  which  0.25  is  the  result  of  combining  several  constants.  Now,  the 
so-called  Corrsin  relation  is  obtained  from  3,14  and  3.16  (Tatarskii 
1971)  . 

C2  -  4  /?T  xx  t'1/3  (3.17) 


P  has  been  determined  for  temperature,  . 8 ,  and  velocity. 


/?u-0.5  (Champagne  et  al .  1977;  Panofsky  and  Dutton  1984).  £  is 

assumed  to  be  equal  to  /9_ .  Since  xu”€  for  velocity,  a  simplified 
Corrsin  relation  is  obtained 


4  p  e' 
'u 


(3.18) 


Using  this  relation,  it  is  simple  to  calculate  e,  the  turbulent 
kinetic  energy  dissipation  rate. 

€  is  an  important  term  in  the  variance  budget  equations  for  a 
turbulent  layer.  These  equations,  neglecting  transport  (Ottersten 
1969),  are  (Fairall  and  Markson  1984) 


de/dt  +  u'w'  du/dz  -  g/6  w —  -e 


2  da^/dt  +  w' T '  de/dz  -  -xr  -  -C*  £1/3/  (4  0T) 


(3.19) 


(3.20) 


where  primed  quantities  are  fluctuations,  barred  quantities  are  means, 
2  2  2 

e-l/2(&u  +ctv  +crw  )  is  the  turbulent  kinetic  energy,  g  is 
the  acceleration  due  to  gravity,  z  is  altitude,  a ^  is  variance  of 
potential  temperature,  Xp  is  the  rate  of  dissipation  of  one  half  the 
thermal  variance,  and  u,  v,  and  w  are  the  horizontal  and  vertical 
velocity  components  respectively.  If  the  flux  Richardson  number  is 
defined  as 


lf  =  [(w '6'  g/6  )/(  u'w'  du/dz ) ] 


and  the  eddy  diffusion  coefficient  is  invoked 


w -  -Kh  86/3 z 


u'w'  -  -K  du/dz 
m 


(3.21) 


(3.22) 


(3.23) 


S 

I 


I 


I 


then  the  following  ratio  results  when  using  the  Corrsin  relations  in 
a  ratio  of  equation  3.20  to  3.19.  (Fairall  and  Markson  1984) 

Cj  /  C2  -1.6  (fl/g)  39/dz  Rf/(1  -  Rf )  (3.24) 

C2  /  C2  -  1.6  (6/ g)  39/dz  Rt/(Pr-  Rt)  (3.25) 

where  R^.  is  replaced  by  the  gradient  Richardson  number, 

Ri  -  g  dlnB/dz  /  (du/dz)2  (3.26) 

and  the  turbulent  Prandtl  number, 

Pr“  V*H  (3-27) 

In  an  actively  turbulent  layer  in  the  free  troposphere , 

R^  «  R^cr^t”0.25  and  P^^l  (Dewan  1980),  then 

C2  /  C2  -  1.6  (9/ g)  89/d z  Rl  (3.28) 

Gradient  Richardson  number  is  the  basis  for  a  model  developed  by 

2 

VanZandt  et  al .  (1978,  1981).  This  model  predicts  -optical  and 

2 

e  based  on  measured  values  of  N  (=g /6  39/dz )  and  S  (=  3u/dz  + 

2 

3v/3z  ) ,  and  input  values  of  L,  the  outer  scale  or  overturning  size 

of  the  turbulent  eddies.  Note  that  the  gradient  Richardson  number  is 

equal  to  N/S.  Based  on  equation  3.28,  this  model  can  also  be  related 

2  2 

to  the  ratio  of  C_  and  C  The  model  is  based  on 

T  u 

Tatarskii's  dimensional  relationships 


16 


where  a  and  b  are  constants.  A  statistical  integration  is  performed 
assuming  that  turbulence  occurs  in  discreet  layers  with  R^<  0.25. 
Theoretical  probability  density  functions  are  assumed  for  shear  and 


temperature  gradients.  The  internal  model  variable  L  is  adjusted  to 


give  best  match  to  radar-derived  profiles.  The  model  is 


inherently  based  on  averaging  over  a  length  much  greater  than  the 
turbulent  scale,  of  which  only  a  fraction,  F,  is  actively  turbulent 
(VanZandt  et  al .  1981).  Thus  high  resolution  in  situ  measurements 
(such  as  measurements  by  aircraft)  may  differ  from  model  predictions. 


3.4  The  Turbulent  Kinetic  Energy  Dissipation  Rate 


Few  in  situ  measurements  of  C  or  e  have  been  made  above  the 
-  u 


2  2 

PBL.  Barat  and  Bertin  (1984)  measured  and  with 


stratospheric  balloons,  and  Fairall  and  Markson  (1984)  with  an 
aircraft.  Kennedy  and  Shapiro  (1980)  calculated  various  turbulence 
parameters  from  research  aircraft  observations  of  CAT  episodes 


associated  with  jetstreams. 


Substituting  for  C  in  equation  3.22  gives  another 


relationship  for  e 


=  3.2  (ff/g)  dd/d z  Ri  e 


2/3 


(3.31) 


This  can  be  used  to  calculate  e  from  C^,  measurements,  and  from 


high  altitude  (where  water  vapor  contribution  is  negligible)  profiler 


measurements  of  -optical  (Tatarskii  1971)  with 


; 79  X  10'6  P/T2)]2 


3.2  (0/g)  dd/dz  Ri  e 


2/3 


(3.32) 


pononur*. 


vt  vx  '.Tremvur.^mvj sri^Tvt^n^n  'vvv'”  ■• 


or  by  assuming  a  Corrsin  relation  for  C  ^ 


C2  -  3.2  *n  £  1/3 


and  assuming 


v  =  K  (dn/dz)‘ 
n  n  ' 


KH  -  R.  N  , 


K  -  K„ 
n  H 


M2  -  (3n/3z) 2 


.  *r^  y-„"»v  y,  'rj  vv  ^.vji 

* 
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(3.33) 


(3.34) 

(3.35) 

(3.36) 

(3.37) 


I 

tr 


i 


then 

C2  -  3.2  R-  M2  N"1  «2/3  (3.38) 

n  l 

where  N  is  Brunt-Vaisala  frequency  squared  and  M  is  Tatarskii's 
refractivity  gradient  (Barat  and  Bertin  1984,  Panofsky  and  Dutton 
1984).  Notice  that  equation  3.28  implies  the  factor  R^  should 
actually  be  R^/(?r  -  R^) ■ 

3.5  Turbulence  Size  Scales 


For  oceanic  turbulence,  Gregg  (1987)  defines  a  turbulence 

3/2  1/2 

bandwidth  based  on  the  ratio  of  Ozmidov  scale  Lq=[€/N  '  ]  '  (a 

buoyancy  scale  at  which  buoyancy  equals  inertial  force)  and  the 

3  1  /  ^ 

Kolmogorov  scale,  r}=[u/e  j  (the  inner  scale  or  cutoff 
wavelength  for  the  inertial  subrange  where  viscous  equal  inertial 
forces) . 


L Q/r,  -  [£AN]3/4 


► 


► 


(3.39) 
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where  1/  is  the  kinematic  viscosity.  The  eddy  overturning  scale  is 
proportional  to  the  buoyancy  scale  in  actively  turbulent  layers. 

Thus,  the  inertial  subrange  occurs  for  size  scales  less  than  the 
buoyancy  scale  and  greater  than  the  inner  scale.  Gregg  provides  the 
following  empirical  guideline  to  interpret  the  bandwidth  ratio 
parameter  as  a  turbulence  activity  parameter. 

Value  of  e/VN  Turbulence  State 

<15  Decaying  turbulence,  fluxes  negligible 

>200  Fully  isotropic 

>10000  Fully  developed  flow 

A  spectrum  begins  to  deviate  from  -5/3  slope  at  size  scales  one  order 
of  magnitude  greater  that  the  Kolmogorov  microscale  (Hill  and  Clifford 
1978).  For  Lo=10rj ,  no  inertial  subrange  should  be  apparent. 
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CHAPTER  4 


THE  EWAK  EXPERIMENT 


4 . 1  Overview 


Two  segments  of  the  EWAK  experiment  were  analyzed  in  this  thesis 
study.  The  first  and  major  segment  was  performed  from  14  April  1986 
to  6  May  1986.  This  segment  was  conducted  at  the  Meteorological  Field 
Site,  located  several  miles  southwest  of  State  College,  PA,  on  the 
Rock  Springs  Agronomy  Research  Farm.  The  relevant  instruments 
operating  at  various  times  during  the  experiment  were  an  Air  Force 
Geophysical  Lab  (AFGL)  scintillometer,  Rome  Air  Development  Center 
(RADC)  scintillometer,  an  AFGL  thermosonde,  an  Air  Force  Weapons  Lab 
(AFWL)  isoplanometer ,  a  Naval  Postgraduate  School  (NPS)  rQ 
scintillometer,  a  Penn  State  University  (PSU)  Doppler  radar  profiler, 
a  PSU  Doppler  sodar,  and  the  Airborne  Research  Associates  (ARA) 
aircraft.  A  summary  of  the  database  relevant  to  this  study  is  given 
in  Table  1. 

The  aircraft  was  based  at  University  Park  Airport,  several  miles 
northwest  of  State  College,  PA.  Data  gathering  flights  were  made  over 
the  field  site  and  insofar  as  possible  within  a  radius  of  10  miles. 
Eight  flights  were  made,  six  yielded  useable  results:  flight  1  was  a 
test  flight  and  flight  5  had  a  total  failure  of  the  instrument 
recorder.  Each  flight  consisted  of  a  fairly  rapid  climb  to  maximum 
altitude  (approximately  30  kft  or  10  km)  and  then  a  slow  descent  (500 
fpm) .  All  descents  were  analyzed  (because  of  the  greater  resolution) . 
but  for  comparison  purposes  several  ascents  were  also  analyzed.  Level 
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Scintillimeter 

Aircraft  GMT  Radar  Wind  Thermosonde  Run  Time  Radar 
Flight  Site  Date  Type  Time  Observ  Time  Flight  Time  AFGL  RADC  Obser 
1  PSU  4/13/86  test  flight 
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flight  data  of  return  trips  by  the  aircraft  to  Boston  was  also  added 
to  several  flights. 

The  second  segment  of  EWAK  was  conducted  at  Griffiss  Air  Force 
Base,  Rome,  NY,  from  11  August  through  15  August  1986.  Both  ARA 
aircraft  and  RADC  scintillometer  measurements  were  taken,  however, 
scintillometer  data  taken  during  these  flights  was  unavailable  for 
this  thesis.  One  additional  flight  was  made  on  20  October  1986.  For 
these  flights,  the  aircraft  was  based  at  Oneida  County  Airport  and 
performed  similar  flight  patterns  centered  around  RADC's  field  site  in 
Verona,  NY. 

4.2  Aircraft  Instrumentation 

A  summary  of  the  instrumentation  available  on  the  aircraft  is 
given  by  Fairall  and  Markson  (1984;  1987).  Thus  only  the  relevant 
turbulence  measurements  are  discussed  here.  Temperature  and  velocity 
signals  were  amplified,  filtered,  and  recorded  on  FM  tape.  Recorder 
gain  settings  were  correlated  to  tape  counters  and  recorded  manually. 
Filters  were  set  for  a  bandpass  between  2  and  200  Hz  for  the 
temperature  signal,  and  5  and  200  Hz  for  the  velocity  signal.  The  FM 
tape  recorder  was  set  at  the  15/16  ips  speed,  which  provided  an 
additional  low  pass  filtering  at  256  Hz. 

Temperature  fluctuations  were  measured  with  a  fast  response 
resistance  bridge  driving  microthermal  sensor  probes  as  the 
temperature  sensitive  resistance  elements.  The  bridge  used  was  a 
Thermo  Systems,  Inc.  (TSI)  type  1044,  dc  Wheatstone  bridge  with  a 
frequency  response  of  800  hz .  The  microthermal  probes  were  TSI  type 
1210  with  T1 . 5  configuration  (W  wire,  4.5nm  dia.,  ice  point  resistance 


R^~5fl) .  A  wingtip  boom  carried  matched  sensors  that  were  connected 
to  the  two  bridge  arms.  One  sensor  was  covered  so  that  bridge  signal 
output  was  proportional  to  temperature  fluctuations  sensed  by  the 
exposed  probe. 

Velocity  fluctuations  were  measured  with  hot-wire  anemometry  in  a 
constant  temperature  circuit.  The  bridge  used  was  a  TSI  model  1050. 
Probes  used  were  the  same  TSI  type  1210  T1 . 5  tungsten  sensors  with  a 
50%  overheat. 

Other  measurements  taken  that  were  used  in  this  analysis  were  air 
temperature  (measured  with  a  standard  Rosemont  sensor  and  bridge 
circuit),  dew-point  temperature  (measured  with  a  Cambridge  Systems, 
Inc.  aircraft  chilled  mirror  instrument),  pressure  altitude  (measured 
with  a  Rosemont  sensor) ,  and  radar  altitude  (measured  with  a  Bonzair 
Mark  10  altimeter) . 

4.3  Additional  Instrumentation 

Both  the  RADC  and  the  AFGL  scintillometers  are  instruments  that 
measure  stellar  scintillation  amplitude  fluctuations  based  on  the  NOAA 
model  II  instrument  described  in  detail  by  Ochs  et  al .  (1977). 
Additional  modifications  were  made  to  the  AFGL  scintillometer  as 
described  by  Murphy  and  Battles  (1986) .  The  AFWL  isoplanometer  and 
NPS  rQ  scintillometer  are  instruments  also  based  on  measuring 
stellar  scintillations  and  may  be  regarded  to  be  derivatives  of  the 
NOAA  instrument.  These  instruments  are  described  in  some  detail  by 
Eaton  et  al.  (1985),  and  Stevens  (1985),  respectively. 


The  AFGL  thermosonde  is  essentially  a  special  temperature 

turbulence  measuring  package  attached  to  a  radiosonde.  The  instrument 

package  is  described  in  detail  by  Brown  et  al.  (1982). 

Penn  State's  Doppler  radar  profiler  is,  basically,  of  the  type 

described  by  Strauch  et  al.  (1986).  It  is  a  VHF  (6  m  wavelength) 

pulsed  Doppler  radar  with  a  phased  array  antenna  located  near  McAlevvs 

Fort,  PA,  on  Shantytown  Road,  approximately  10  miles  south  of  State 

College,  PA.  With  the  profiler,  vertical  profiles  of  wind  speed  and 

direction  are  obtained  as  well  as  radar  reflectivity  (backscatter 

cross  section) .  Radar  volume  reflectivity  can  be  used  to  calculate 
2 

Cn  values  at  each  of  the  radar  gates  (VanZandt  et  al .  1978).  The 
profiler  produces  both  high  and  low  resolution  profiles.  High 
resolution  is  -290  m  from  -1  km  to  8  km  MSL.  Low  resolution  is  -870  in 
from  -1.5  km  to  17.5  km  MSL. 

PSU  Doppler  sodar  data  was  not  used  in  this  study. 
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CHAPTER  5 


DATA  PROCESSING 


The  relation  of  the  structure  constant  to  the  variance  spectral 

density  as  a  function  of  wavenumber  (Equation  3.16)  allows 

2  2 

determination  of  C^,  and  C^  from  spectral  analysis  of 
turbulence  data.  The  aircraft  data  was  a  time  series  record  of 
voltage  fluctuations,  which  could  be  transformed  into  frequency 
spectra.  Taylor's  frozen  turbulence  hypothesis  allows  for  translation 
of  the  spatial  spectrum  0(k^)  to  a  frequency  spectrum  <£'(f) 

(Panofsky  and  Dutton  1984) ,  by  the  relation 


k-^  -  27rf/u 


(5.1) 


where  k^  is  the  spatial  wavenumber,  u  is  the  relative  speed  of  the 
passing  eddies  (in  this  case,  u-60  m/s,  the  aircraft  speed),  and  f  is 
the  frequency  in  Hz.  Since  total  spectral  intensity  must  be  invariant 
(Panofsky  and  Dutton  1984) ,  temporal  and  spatial  spectra  can  be 
related  by 


0(k1)dk1  =  4>'  (f)df 


Substitution  gives 


(5.2) 


<^>'(f)  =  0(2?rf/u)dk^  /df  =  <^(27rf/u)  2r/u 


(5.3) 


and  then 


<t>j(  f)  -  2tt/u  <^T(k1) 


^(f)  -  2tt/u  (0.25  C^  k*5/3  ) 


(5.4) 


(5.5) 


^(f)  -  (2*/u)'i/J  (0.25 


(5.6) 


Similarly  for  velocity 

^(f)  -  (2i/u)'2/3  (0.25  C2  f‘5/3) 


(5.7) 


Rearranging  gives  the  equation(s)  relating  the  structure 
functions  to  the  measured  spectra 


CT,u  =  4  (2w/u) 2/3  ^u(f)  f5/3 


(5.8) 


where  <f>' ^  u(f)  is  a  representative  value  in  the  inertial  subrange 
(Fairall  and  Markson  1984) . 

Another  method  to  determine  the  structure  constants  (the  single 
probe  RMS  method)  is  suggested  by  the  relation  of  total  variance  to 
spectral  density  (Fairall  and  Markson  1984) 


Variance  -  mean  square  - 


*'(f)  df 


(5.9) 


where  f^-upper  frequency  limit,  and  f^-lower  frequency  limit. 

Thus  the  filtered  RMS  of  the  signal  can  be  used  to  directly  obtain  th 
structure  constant  by  integrating  5.6  or  5.7 


Variance  =  (RMS)2=  (2ir/u) '2/30. 25  C2  u(-  \  )(fu’2/3'  f^2^) 

(5.10) 


Rearranging  gives 


r2  8  ,0  ,  ,2/3  ,DMC , 2  ..-2/3  .-2/3  ,-l 

CT  u  ”  3  (2tt/u)  (RMS)  '  fu  ) 


(5.11) 


2 


The  thermosonde  instrument  calculated  directly  using  the 


double  probe  RMS  method.  If  two  probes  with  a  known  separation  /  are 
used,  then  with  an  RMS  of  the  difference  of  the  signal 


(Brown  et  al.  1982) 


CT  ”  <  (T1  '  T2)2  >  /  ^/3  "  (^S)2/  ?2/3  (5.12) 

The  spectral  method  has  the  advantage  of  greater  accuracy  over 
the  RMS  method  because  a  region  of  the  spectrum  can  be  selected  which 
is  relatively  free  of  noise.  The  RMS  variance  includes  all  noise  in 
the  bandwidth  analyzed  (Fairall  and  Markson  1984).  However,  the 
disadvantage  to  the  spectral  method  is  the  greatly  increased 
computation  time.  RMS  methods  can  be  done  on  a  virtually  real  time 
basis . 

The  aircraft  temperature  turbulence  data  was  taken  as  voltage 
fluctuations  in  bridge  output  due  to  resistance  fluctuations  of  the 
probes.  Temperature  fluctuations  cause  these  resistance  changes. 

This  relation  is  given  by 

dV/dt  -  i  Gfi  dR/dt  (5.13) 

where  the  current  to  the  probe  sensor  is  i,  and  bridge  amplification 
or  gain  is  Gg(=4000)  (Fairall  and  Markson  1984) . 

The  resistance-temperature  relation  is  linear  for  tungsten 
sensors  over  a  broad  temperature  range ,  giving 

R  =  R.[l  +  a(T  -  T.)]  (5.14) 

where  R^  and  T^  are  values  at  the  ice  point  (Ri=5Q) ,  and  a  is 
the  effective  thermal  resistance  coefficient  (a=.0026  K  ^)  ,  which 
includes  the  loss  of  sensitivity  due  to  probe  support  effects  (Fairall 
and  Markson  1984).  Thus, 

dR/dt  -  R,.  a  dT/dt 


(5.15) 


(5.16) 


dV/dt  -  i  Gfi  G  (R.  a  dT/dt) 

where  the  added  G  term  accounts  for  any  additional  signal  conditioning 
or  gains  applied  (usually  G=10) .  The  voltage  variance  spectrum  is  the 
Fourier  transform  of  the  square  of  voltage  fluctuations  which  then  can 
be  related  to  the  temperature  variance  spectrum  (Fairall  and  Markson 
1984) , 

^  -  (i  G  Gg  R.  a)2  ^(f)  (5.17) 

Including  noise  sources,  this  would  be 

4>lj  -  (i  G  Gfi  Rt  a) 2  +  Nu  +  Nt  +  Nv/(iR.a)2  )  (5.18) 

where  is  broadband  voltage  noise  in  the  system,  NT  is 
temperature  noise  detected  by  the  sensors  but  not  due  to  turbulence 
(e.g.,  sonic  adiabatic  compression),  and  is  velocity 
contamination  due  to  cooling  rate  variations  associated  with  velocity 
fluctuations.  An  indepth  study  of  signal-to-noise  ratio  was  done  by 
Fairall  and  Markson  (1984).  It  was  concluded  that,  in  flight,  most 
noise  was  temperature  noise.  Although  the  source  of  this  noise  could 
not  be  identified,  it  could  be  lessened  by  changing  sensor  location. 
Optimal  sensor  location  proved  to  be  on  the  wingtip.  Sensor  currents 
were  kept  at  a  maximum  just  below  the  velocity  contamination 
threshhold  (i=5  mA) . 

Velocity  turbulence  data  was  also  taken  as  a  time  series  of 
voltage  fluctuations.  King's  Law  voltage  relationship  (Hinze  1975) 


2 

where  V  and  B  are  calibration  constants  (Fairall  and  Schacher 
o  v 

2  2 

1977)  computed  for  50%  overheat  at  STP  as  VQ  =2.7  Volts  and 
2  -11/2 

Bq=1 . 0  Volts  /(ms  )  '  .  Correcting  for  varying  temperature 

and  pressure  gives 

B  -  Bq  (P/1013)1/2  (288/T)1/4  [ (453 -T)/288 ]  (5.20) 

(Fairall  and  Markson  1984) .  Velocity  fluctuations  cause  bridge  outpu 
voltage  variations  as 

dV/dt-  [  B  G/(4  V  u1/2)]  du/dt  (5.21) 

with  G  additional  signal  conditioning  gains  (usually  G-100)  and  V  the 
mean  bridge  voltage  (about  3  Volts  for  nominal  aircraft  speeds) . 

Relating  the  square  of  fluctuations  to  the  variance  spectrum 
gives 

<t>y  -  (B  G/  4  V  Vu)2  ^  (5.22) 

2  2 

The  final  forms  used  for  calculating  and  from  the 

aircraft  data  (Fairall  and  Markson  1984)  were 

C2  =  4  ( 27r/u) 2/3  (i  G  Gg  R.  a)'2  ^(T)(f)  fV3  (5.23) 

C2  -  4  ( 2tt ) 2/3  UV3  (4  v/  B  g)2  ^(u)(f)  f5/3  (5.24) 

Flight  12  gain  settings  were  not  recorded.  They  were  assumed  to 
be  the  same  as  flight  11  for  processing,  but  for  that  flight,  this 


should  be  noted. 


Initially,  the  temperature  and  velocity  channels  of  the  taped 
data  were  converted  from  analog  to  digital  with  an  Infotek  A/D  board. 
Approximately  two  minutes  of  analog  data  were  digitized.  A  128  point 
Fast  Fourier  Transform  (FFT)  was  done  and  a  final  spectrum  produced 
from  an  average  of  400  spectral  blocks.  This  gave  a  vertical 

resolution  of  0.3  to  0.5  km,  about  the  same  as  the  radar  profile. 

2 

The  log  of  power  spectral  density  (Volts  /Hz)  was  plotted  as 
2 

decibels  [10*log  (V  /Hz)]  against  log  of  frequency  (Hz)  for  4  to  256 
Hz.  A  key  frequency  was  chosen  that  consistently  appeared  in  the 
inertial  subrange  (in  a  region  of  high  signal  to  noise)  and  the 
spectral  density  at  this  point  was  used  to  plot  a  -5/3  slope  through 
the  spectrum  to  show  goodness -of-f it  to  the  Kolmogorov  turbulence 
theory.  Due  to  the  small  number  of  points,  the  spectra  obtained  with 
this  method  were  somewhat  crude  and  very  difficult  to  interpret. 

Figure  1  shows  three  different  spectra  from  flight  2  as  examples 
(indicated  as  EWAK2BT-rec#) .  Record  1  has  a  high  signal-to-noise 
ratio,  record  31  is  very  noisy  with  little  signal,  and  record  45  is  an 
unusual  steeply  sloped  spectra. 

A  second  method  gave  far  more  satisfactory  results.  Roughly  two 
minutes  of  data  were  analyzed  by  a  dual  channel  Hewlett-Packard  3562A 
dynamic  signal  analyzer  by  linearly  averaging  55  times  with  an  overlap 
of  75%.  A  1600  point  FFT,  using  a  Hanning  window,  produced  a  power 
spectrum  from  .125  to  100  Hz.  This  frequency  range  corresponds  to  a 
wavenumber  range  of  0.01  to  10  m  ^ .  Resolution  of  the  spectra  was 
greatly  increased.  The  same  three  records  from  flight  2  are  shown  in 
Figure  2. 
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Figure  1.  Voltage  variance  spectra  from  first  FFT  program  (128 
points) .  For  each  record,  the  top  spectrum  (channel  wl)  is 
temperature  variance  and  the  bottom  spectrum  (channel  ~2 )  is  velocity 
variance.  Voltage  power  spectral  density  (dB)  is  plotted  against 
frequency  (log  Hz).  The  slanted  line  represents  a  -5/3  slope  for 
comparison.  Coded  flight  name  (e.g.,  EWAK2BT)  and  record  number 
(e.g.,  -  1),  tape  footage,  and  altitude  are  shown  interior  to  each 
channel  #1  graph.  Three  different  dual -channel  records  in  flight  2 
are  shown’  1,  31,  and  45. 
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Figure  2.  Voltage  variance  spectra  produced  by  the  signal  analyzer's 
1600  point  FFT.  For  each  record,  Ovl  (top)  is  temperature  variance 
and  Ov2  (bottom)  is  velocity  variance.  Voltage  power  spectral  density 
(dB)  is  plotted  against  frequency  (log  Hz) .  Processing  parameters  are 
stated  above  each  spectrum.  Coordinates  of  the  cursor  are  given  in 
the  upper  left  corner  of  each  spectrum.  Vertical  axis  range  has  been 
fixed  as  stated  in  the  lower  left  of  each  spectrum.  Coded  flight  name 
(e.g.,  EAIR)  and  record  number  (e.g.,  -  1),  tape  footage,  and  altitude 
are  annotated  above  each  Ovl  graph.  The  same  three  dual-channel 
records  in  flight  2  are  shown  as  in  Figure  1. 
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CHAPTER  6 


ANALYSIS  METHODS 

6.1  Aircraft  Turbulence  Data  Analysis 

An  indepth  analysis  on  aircraft  instrumentation/experimental 
error  was  done  by  Fairall  and  Markson  (1984) .  Uncertainty  due  to 
measurement  error  will  therefore  not  be  discussed  here.  However,  for 
this  thesis,  much  effort  was  put  into  detailed  analysis  of  the 
aircraft  turbulence  data  in  attempts  to  improve  analysis  precision. 
These  processes  are  described  in  this  chapter. 

To  quantify  goodness-of-fit,  a  linear  regression  analysis  was 
performed  on  the  log- log  spectral  data  and  percent  error  was 
calculated  for  deviation  from  -5/3  slope.  Most  spectra  exhibited  the 
-5/3  dependence  out  to  the  broadband  noise  level  of  the  instruments, 

»3  X  10  ^  K^m  for  temperature  and  =1  X  10  ^  (m/s)^m  for 
velocity  (Fairall  and  Markson  1985).  Because  of  this  background  noise 
level,  notable  in  the  high  frequency  range  (over  20  Hz),  it  was 
determined  that  the  best  fit  would  be  obtained  at  the  lowest  frequency 
possible  within  the  inertial  subrange.  Unfortunately,  the  temperature 
data  had  been  low  pass  filtered  at  2  Hz  and  the  velocity  data  at  5  Hz. 
Filters  of  the  type  used  on  the  aircraft  (Thermo  Systems,  Inc.  model 
1057  signal  conditioner  with  TSI  model  1051-2  monitor  and  power 
supply)  were  analyzed  with  white  noise  input  on  the  signal  analyzer. 
Then  a  power  spectrum  ratio  of  unfiltered  noise  to  filtered  noise  was 
calculated  (Figures  3a  and  3b).  This  function  could  then  be 
multiplied  point  by  point  to  essentially  "defilter"  the  data. 


FILTER  FUNCTION  VALUE 


Spectral  drop-off  was  still  noted  in  the  low  frequency  end  (<1 
Hz)  of  the  spectra  after  this  filter  function  was  applied.  Several 
reasons  for  this  drop-off  were  proposed:  a  real  phenomenon  (i.e.,  the 
outer  scale  of  the  inertial  subrange  or  the  production  size  of  the 
turbulence),  a  non-linear  characteristic  of  the  filter  when  the  signal 
is  log-distributed  and  not  linearly  distributed  like  the  white  noise, 
or  individual  error  in  the  filter  (two  identical  filters  were  tested 
and  found  to  have  different  signatures) .  Assuming  f  is  1  Hz  and  u  is 
60  m/s,  the  frozen  turbulence  theory  (A  ^-27rf/u)  gives  A  of  10  m  as 
the  production  scale.  For  the  altitudes  flown,  this  is  within  an 
order  of  magnitude  of  typical  estimates  (VanZandt  et  al .  1981;  Gregg 
1987;  Barat  and  Bertin  1984).  For  comparison  to  the  characteristics 
of  the  measured  filter  function,  an  atmospheric  filter  function  was 
derived  to  compensate  for  non-linear  filter  effects  and  any  unknown 
atmospheric  effects  (e.g.,  variable  production  size  scale).  Several 
sample  spectra  from  the  boundary  layer  with  exceptionally  high 
signal- to-noise  ratio  and  -5/3  slope  were  used  to  derive  a  function 
based  on  -5/3  slope  (Figures  4a  and  4b). 

All  spectra  were  then  run  through  a  three  pass  progressive 
editing  program.  If  the  percent  error  on  a  regression  fit  from  5  to 
18  Hz  was  greater  than  20%,  the  spectrum  was  rejected  on  pass  1. 

Figure  5a  again  shows  records  1,  31,  and  45  of  flight  2  for  pass  1. 

On  the  second  pass,  the  atmospheric  filter  was  used  to  defilter  the 
data  from  0.1  to  5.6  Hz  (the  cutoff  frequency  where  both  2  and  5  Hz 
filters  were  approximately  equal  to  one).  Figure  5b  shows  the  same 
three  records  for  pass  2.  If  percent  error  was  greater  than  15%  on  a 
second  regression  fit  from  2  to  18  Hz,  then  the  spectrum  was  run 
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Figure  5a.  Voltage  variance  spectra  showing  editing  procedure.  These 
spectra  are  filtered  spectra  derived  from  the  taped  data.  This 
represents  editing  pass  1.  Voltage  power  spectral  density  (dB)  is 
plotted  against  frequency  (log  Hz).  For  each  record,  channel  1  (top) 
is  temperature  variance  and  channel  2  (bottom)  is  velocity  variance. 
Channel  number,  pass  number,  and  coded  flight  name  and  record  number 
are  shown  interior  to  each  graph.  The  slanted  line  represents  the 
regression  fit.  Regression  slope,  correlation  coefficient,  and 
percent  error  are  given  in  the  upper  left  corner  of  each  spectrum. 

The  same  three  dual-channel  records  in  flight  2  are  shown  as  in  Figure  1. 
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Figure  5a.  (continued) 


Figure  5a.  (continued) 
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Figure  5b.  Voltage  variance  spectra  showing  editing  procedure.  The 
spectra  have  been  defiltered  with  the  atmospheric  filter  function. 
This  represents  editing  pass  2.  Voltage  power  spectral  density  (dB) 
is  plotted  against  frequency  (log  Hz).  For  each  record,  channel  1 
(top)  is  temperature  variance  and  channel  2  (bottom)  is  velocity 
variance.  Channel  number,  pass  number,  and  coded  flight  name  and 
record  number  are  shown  interior  to  each  graph.  The  slanted  line 
represents  the  regression  fit.  Regression  slope,  correlation 
coefficient,  and  percent  error  are  given  in  the  upper  left  corner  of 
each  spectrum.  The  same  three  dual-channel  records  in  flight  2  are 
shown  as  in  Figure  1. 


Figure  5b.  (continued) 
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Figure  5c.  Voltage  variance  spectra  showing  editing  procedure. 
Average  noise  level  has  been  subtracted  from  each  defiltered  spectrum 
This  represents  editing  pass  3.  Voltage  power  spectral  density  (dB) 
is  plotted  against  frequency  (log  Hz).  For  each  record,  channel  1 
(top)  is  temperature  variance  and  channel  2  (bottom)  is  velocity 
variance.  Channel  number,  pass  number,  and  coded  flight  name  and 
record  number  are  shown  interior  to  each  graph.  The  slanted  line 
represents  the  regression  fit.  Regression  slope,  correlation 
coefficient,  and  percent  error  are  given  in  the  upper  left  corner  of 
each  spectrum.  The  same  three  dual-channel  records  in  flight  2  are 
shown  as  in  Figure  1. 
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through  pass  3.  An  average  noise  level  was  calculated  from  the  twenty 

lowest  spectral  density  values  and  subtracted  from  the  entire 

spectrum.  Figure  5c  shows  the  same  three  records  for  pass  3.  If  the 

percent  error  on  a  new  regression  fit  from  2  to  18  Hz  was  still 

greater  than  15%,  the  spectrum  was  rejected  on  pass  3  and  flagged. 

Spectra  accepted  on  one  of  the  above  three  passes  were  coded  as  "1". 

Flagged  spectra  were  then  manually  examined.  If  a  spectrum  contained 

a  visible  portion  of  the  inertial  subrange  that  was  represented  by  the 

regression  fit,  it  was  coded  "2".  Spectra  with  too  small  a 

signal- to-noise  ratio  (no  visible  inertial  subrange)  were  coded  "3". 

Unreal  spectra  were  coded  as  "4".  These  unreal  spectra  were  steeply 

sloped  in  log- log  space.  No  theory  predicts  such  a  spectrum.  Thus, 

these  spectra  were  thought  to  be  byproducts  of  measurement  error, 

either  overdriving  the  filter  (i.e.,  amplitude  fluctuations  so  great 

that  filter  saturation  occurred)  or  noise  distortion  (e.g.,  change  of 

aircraft  power  setting  or  a  radio  transmission) .  Code  4  spectra  were 

eliminated  from  the  analysis.  All  other  analyses  was  performed  on 

code  1  and  2  spectra,  unless  otherwise  noted.  Table  2  gives  a 

representative  breakdown  of  spectra  by  editing  codes. 

Based  on  the  regression  analysis,  the  following  formulas  were 

used  to  derive  the  spectral  density  necessary  to  calculate  the 

2  2 

structure  parameters,  and  .  'X'  was  assigned  to 

2 

frequency  and  'Y'  to  Volts  /Hz  (spectral  density).  Using  the 
power-law  equation, 


Y  -  a  Xm 
x5/3 


(6.1) 


Y 


yOn  +  5/3) 


/  A  9  ^ 


U,' 


Table  2. 

Representative  editing  results  for  each  flight  and  in  total.  The 
breakdown  of  spectra  in  each  pass  for  code  1  is  shown,  as  well  as  the 
number  of  spectra  in  the  remaining  codes. 
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Table  2. 
(continued) 


Velocity 

Flight 

Spectra 
Code  1  : 
Pass  1 

2 

3 

Code  2 

2 

27 

4 

9 

2 

3 

21 

0 

4 

1 

4 

6 

0 

0 

0 

6 

34 

0 

3 

7 

71 

25 

0 

0 

1 
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36 

0 
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0 

8 

28 

2 

3 

2 

9 

28 

0 

3 

1 

10 

53 

1 

0 

0 

111 

33 

0 

2 

2 

1111 

16 

0 

2 

2 

12 

21 

1 

5 

2 

TOTALS 

Temperature  Spectra 

Code  1  :  Pass  1 

2 

3 

Code  2 

222 

25 

55 

1 

Velocity 

Spectra 

Code  1  : 

Pass  1 

2 

3 

Code  2 

328 

8 

33 

2 

q  -  x5/3  y 

a  -  Y  /  Xm 

Q  -  10<<1oS  '  m  <lo§  x>>  i0(m  +  5/3)  <loS  *> 

Q  _  1Q(<log  Y>  +  5/3  <log  X» 

Q  _  1Ql(S  log  Y)/N  +  1.667  (S  log  X)/N) 


(6.3) 

(6.4) 

(6.5) 

(6.6) 

(6.7) 


where  Q  is  the  spectral  density  times  frequency  to  the  5/3  and  is 

based  on  the  spectral  density  and  frequency  values  midrange  in  the 

regression  fit.  Q  is  a  necessary  input  parameter  for  the  C^'  and 

2 

equations  (equations  5.23  and  5.24),  calculated  from  the 
temperature  and  velocity  turbulence  spectra,  respectively. 

6.2  Mean  Aircraft  Meteorological  Profiles 

The  aircraft  temperature  and  humidity  profiles  were  digitized  at 
levels  corresponding  to  the  profiler  high  resolution  range  gate 
heights.  Dewpoint  temperature  was  converted  to  vapor  pressure,  e, 
over  water  using  the  following  formula  (Iribarne  and  Godson  1981). 

log  e  =  -(2937.4  /  T)  -  4.9283  log  T  +  23.5471  (6.8) 

where  e  is  vapor  pressure  in  mb  and  T  is  in  Kelvin.  Vapor  pressure 
was  then  converted  to  specific  humidity,  q,  using  the  approximation 
i  e  =  (P  q)/0.622.  Potential  temperature  and  temperature  have  the 

defined  relationship 


ff  =  T  (1000  mb/P) 


0.286 


(6.9) 


with  pressure  from  the  standard  atmosphere  relationship 

P  -  P  e'(z/8)  (6.10) 

o 

where  Pq  is  1013.2  mb,  altitude  z  is  in  km  and  scale  height  H  is  =8 
km  (Iribarne  and  Godson  1981;  Wallace  and  Hobbs  1977).  Aircraft 
potential  temperature  profiles  are  shown  for  all  flights  in  Figure  6. 
Aircraft  specific  humidity  profiles  for  only  the  PSU  flights  are  shown 
in  Figure  7. 

6 . 3  Other  Data  Sources 

Profiler  data  was  logged  at  one  hour  intervals  denoted  by  the  GMT 

time  (zulu  time,  "z")  at  the  end  of  the  average  period.  Height  is 

2 

given  in  m,  wind  speed  in  m/s,  direction  in  degrees,  and 
-2/3 

(radar)  in  m  (Mike  Moss,  letter  to  author,  July  1987;  Mike  Moss, 
letter  to  author,  November  1987).  Radar  wind  profiles  for  the  PSU 
flights  are  shown  in  Figure  8. 

AFGL  thermosonde  data  was  available  at  20  m  resolution.  The 

following  information  was  available  for  each  level:  altitude  (km), 

pressure  (mb),  temperature  (°C) ,  relative  humidity  (%) ,  C  ^ 

-2/3 

(m  ),  wind  speed  (m/s),  and  direction  (degrees)  (Robert  Beland, 

letter  to  William  Syrett,  September  1986).  AFGL  thermosonde  potential 

temperature  profiles  corresponding  to  aircraft  flights  are  shown  in 

Figure  9.  Corresponding  AFGL  thermosonde  specific  humidity,  mean 

2 

wind,  and  calculated  raw  profiles  are  shown  in  Figures  10,  11, 

and  12,  respectively. 

AFGL  scintillometer  data  was  provided  as  seven- level  profiles 
2  -2/3 

with  in  m  ,  and  height  (corrected  for  zenith  angle)  in  km 
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Figure  6.  (continued) 
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Figure  8.  (continued) 
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Figure  12.  AFGL  thermosonde-cal^ulated  vertical  profiles  of 
C  -optical  based  on  measured  C^.  .  Altitude  is  MSL. 
tR  ermosonde  flight  number  is  given  interior  to  each  graph. 
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(Murphy  and  Battles  1986) .  AFGL  scintillometer  profiles  corresponding 

to  aircraft  flights  are  shown  in  Figure  13. 

RADC  scintillometer  data  was  provided  as  seven- level  profiles 
2  -2/3 

with  C  in  m  '  at  the  standard  kernel  levels  (AGL  altitude 
n  v 

assuming  zenith  angle  of  zero)  (Donald  Stebbins,  letter  to  author, 
September  1986).  Only  one  data  run  corresponded  to  the  aircraft 
flight  times.  It  is  shown  in  Figure  14. 

NPS  rQ  scintillometer  data  was  provided  as  rQ  in  cm  as  a 
function  of  Zulu  time  (Donald  Walters,  telephone  conversation, 

November  1987) .  Data  taken  during  aircraft  flight  was  limited  to  one 
night,  4  May  1986.  The  average  value  of  rQ  during  this  time 
segment,  0400-0540  Z,  was  5  cm. 

AFWL  isoplanometer  data  was  unavailable  for  use  in  this 


comparison. 


CHAPTER  7 


RESULTS  AND  DISCUSSION 
7.1  Analysis  of  C^/C^ 

2  2 

The  vertical  profiles  of  C^  and  C^,  calculated  from  the 

aircraft  data  show  strong  correlation  for  flights  2,  6(ascent  and 

descent),  8(ascent  and  descent),  9(ascent),  lO(descent) ,  11,  and  12 

(Figure  15) .  The  other  flights  also  show  some  correlation,  but  not  at 

all  levels.  Correlation  of  higher  values  of  these  turbulence 

parameters  indicates  regions  of  the  atmosphere  that  contain  turbulent 

2  2 

layers  (Ottersten  1969) .  C^,  and  C^  are  directly  correlated 

o 

at  all  levels  for  all  flights,  as  is  to  be  expected  since  C^  is  a 

2 

function  of  C^,  . 

A  rough  estimate  of  boundary  layer  height  can  be  determined  by 

spectral  signal-to-noise  ratio  for  these  profiles,  since  the  signal  is 

much  stronger  in  the  turbulent  boundary  layer.  For  this  purpose 

spectra  were  examined  in  detail  manually.  Table  3  gives  an  estimate 

2  2 

of  PBL  height  for  each  flight.  If  C^,  versus  C^  is  graphed 
(for  values  above  the  boundary  layer  where  the  TKE  and  variance  budget 
approximations  discussed  in  Chapter  3  are  valid),  a  power-law  relation 
is  evident  (Figures  16  and  17).  This  ratio  is  proportional  to 
Richardson  number  and  temperature  gradient  (right  hand  side  of 
equation  3.28).  All  flights  except  6(ascent  and  descent),  7 (descent), 
and  ll(ascent)  show  distinct  sloped  correlations.  All  flights  were 
plotted  simultaneously  and  a  regression  analysis  was  done  to  determine 
slope  and  goodness-of - f it  (Figure  18).  The  plot  of  flight  6  had  high 
scatter.  Removal  of  this  flight  decreased  the  scatter  and  thus 
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Figure  15.  (continued) 
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Figure  16.  (continued) 
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increased  correlation  of  the  entire  data  base  (Figure  18) .  The  slope 

for  all  the  vertical  profiles  was  close  to  3/5.  Correlation  and  slope 

for  all  the  level  flights  showed  approximately  the  same  scatter 

(correlation)  and  slope  as  the  vertical  profiles  (Figure  19) . 

2  2 

In  Figure  20,  C^,  was  plotted  versus  Cu  multiplied  by 

9/ g  ( ae/dz )  -  T/g  (dT/dz  +  n  (7.1) 

for  all  vertical  flights  (having  associated  temperature  profiles) . 

Again,  flight  6  had  high  scatter  and  flight  7  showed  little  slope.  A 

combined  plot  of  all  flights  showed  slightly  higher  correlation 

without  flight  6  and  again  roughly  a  slope  of  3/5  (Figure  20) .  This 

ratio  corresponds  to  values  of  1.6  R^/(Pr*  R^)  that  range  from 

o 

0.3  for  high  values  of  C^,  to  as  much  as  10  for  low  values  of 


2  2 

As  an  additional  comparison,  /  (Cu  9/ g  39/3 z)  was 
plotted  against  a  turbulence  activity  parameter  (Gregg  1987) ,  e/i/N 
(Figure  21) .  N  was  calculated  as  g  31n9/3z  from  a  polynomial  fit  to 
the  aircraft  potential  temperature  profile.  The  formula  used  to 
derive  kinematic  viscosity  was  taken  from  the  U.  S.  Standard 
Atmosphere  Supplements .  1966  (Environmental  Science  Services 
Administration  et  al.  1966). 

M  -  (0  T3/2)/(T  +  S)  (7.2) 

where  p.  is  the  dynamic  viscosity,  $  -  1.458E-6  kg  sec  3K  3/2 . 

and  S  -  110.4  Kelvin.  Kinematic  viscosity,  i/ ,  equals  n/p.  Density 

.  3 

was  calculated  as  p  -  exp(-Z/8  km)  where  pQ  is  1.225  kg  m 

2  2 

T  was  taken  from  the  aircraft  profiles.  C_  ,  C  ,  and  N 


jrrjrrjr:  v  vw  v-j*  v 


Cu-'2»(T/'a)«(dT/dZ  +  r)  (K'-a  m-'-a/'3) 


Slope  -  .554 

Correlation  Coefficient  =  .703 


Figure  20.  A  composite  comparison  of  C^,  versus  C 
normalized  by  the  expression  (T/g  (3T/3z  +  T) ] .  T^e  broken  line 
represents  a  regression  analysis.  Slope  and  correlation  coefficient 
are  given  below  each  graph. 
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were  all  calculated  from  the  aircraft  data.  Figure  22  is  a  composite 

of  all  vertical  flights.  This  plot  clearly  shows  that  as  turbulence 

2  2 

activity  goes  down,  the  ratio  of  to  normalized  goes 

2  2  2 

up.  For  activity  values  above  10  ,  the  C^,  /C  (norm)  ratio 
becomes  constant  at  about  0.4. 

7.2  Mean  Gradients  and  Richardson  Number 

Next,  the  Richardson  number  was  evaluated.  Since  VanZandt  et  al . 

(1981)  defined  R^  as  N/S ,  N  versus  S  was  plotted,  thus  showing 

Richardson  number  as  a  slope  (Figure  23) .  N  was  calculated  by  a 

finite  difference  derivative  between  each  level  in  the  digitized 

aircraft  potential  temperature  profiles,  while  S  was  calculated  by 

finite  difference  from  the  profiler  hourly  average  taken  during  the 

cime  the  aircraft  was  collecting  data  (Table  1) .  S  was  calculated  as 
2  2 

(3u/3z)  +  (3v/3z)  .  Figure  24  is  a  composite  of  all  flights 

shown  in  Figure  23.  Most  values  show  Ri>  0.25  above  the  broken 

line,  with  the  shear  term  on  the  order  of  10  ^  and  the  buoyancy  term 
-  5  -4 

ranging  from  10  to  4  X  10  Richardson  number  profiles  for 
each  flight  were  directly  calculated  (Table  4) ,  and  then  averaged 
(Table  5 ) . 

These  values  of  N  and  S  are  the  necessary  input  parameters  for 

the  VanZandt  et  al.  model  which  was  used  to  predict  vertical  profiles 
2 

of  C  and  e .  This  was  done  for  each  flight  conducted  at  Penn 
n  ° 

State  where  profiler  data  was  available  (Figure  25).  The  model  ourr 

2 

is  normalized  optical  .  Standard  pressure  calculated  from  the 
geometric  height  and  the  aircraft  temperature  profile  were  used  to 
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Table  4. 

Calculated  vertical  profiles  of  Richardson  number  based  on  R.— N/S 
and  using  N  and  S  values  shown  in  Figure  17.  Altitude  shownXis  MSL 
and  is  the  average  altitude  between  each  vertical  radar  observation 
level . 
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Table  4. 
(continued) 


Flight  4 
Alt  (km) 

Flight  6 
Alt  (km) 

R. 

1.4 

2.45 

1.4 

-  .03 

1.7 

34.24 

1.7 

142.19 

2 

86.56 

2 

.  78 

2.3 

5.99 

2.3 

21.71 

2.6 

5.87 

2.6 

8.65 

2.9 

1.09 

2.9 

179.84 

3.2 

-  .02 

3.2 

1.6 

3.4 

39.34 

3.4 

3.2 

3.7 

3.23 

3.7 

181.64 

4 

6.27 

4 

7.6 

4.3 

6.21 

4.3 

2.24 

4.6 

10.47 

4.6 

9.6 

4.9 

6.66 

4.9 

27.89 

5.2 

3.32 

5.2 

4.13 

5.5 

5.11 

5 . 5 

1.92 

5.8 

1.84 

5.8 

2.19 

6.1 

9.16 

6.1 

.05 

6 . 4 
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6.4 
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6.6 
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Table  4. 
(continued) 


Flight  7 
Alt  (km) 

Ei 

Flight  8 
Alt  (km) 

R. 

“JL 

1.4 
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1.4 

.02 

1.7 
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1.7 

0 
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2.31 

2.3 

1.81 

2.3 

8.51 

2.6 

1.28 

2.6 

8.05 
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.82 
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3.2 
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3.2 
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3.4 

1.67 

3.4 
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3.7 

32.94 

3.7 
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4 
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4 
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4.3 

.53 
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4.54 
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67 . 55 
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5.5 
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6.4 
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6 . 6 
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6.6 

23.27 

6.9 
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Figure  25.  Comparison  of  the  VanZandt  et  al .  model  vertical  profile 
prediction  of  -optical  (solid  line)  and  the  aircraft 
C_  -derived  vertical  profile  of  C  -optical  (broken  line). 

PSU  flights  are  shown  (labeled  interior  to  each  graph) . 
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Figure  25.  (continued) 


”de- normalize"  the  model-derived  parameter.  Since  only  the  high 
resolution  profiler  range  gates  were  used  to  calculate  S  at  a 
resolution  of  290  m,  the  prediction  was  only  calculated  from  1410  m 
MSL  to  7515  m  MSL.  An  average  MSL  altitude  for  State  College,  PA,  is 
300  m.  Except  for  flight  6(ascent)  the  model  and  aircraft  profiles 
agreed  to  within  an  order  of  magnitude.  This  particular  discrepancy 
in  flight  6  is  thought  to  be  the  result  of  errors  in  gain  records  for 
the  flight. 

7.3  Profiles  of  C  ^ 
n 

2 

The  VanZandt  et  al.  model  can  also  be  used  to  give  radar  . 

2 

The  relationship  between  profiler  measured  C -radar  and 
2 

Cn  -optical  is  a  function  of  temperature  and  humidity  (Tatarskii 

2 

1971).  This  parameter,  previously  referred  to  as  X  ,  has  been 

calculated  using  the  aircraft  specific  humidity  and  temperature  data. 

2  2  2  2 
X  is  the  ratio  of  C  -radar  to  C  -optical.  X  at  each 

n  n 

aircraft  average  observation  level  was  calculated  with  a  polynomial 

fit  to  the  temperature  profile  and  a  finite  difference  derivative  from 

the  humidity  profile.  Smoothing  of  the  vertical  profiles  was  achieved 

by  setting  negative  or  zero  temperature  derivatives  equal  to  0.00001 

K/m  and  positive  specific  humidity  profiles  equal  to  zero.  This 

2 

calculation  is  given  in  Table  6.  X  at  each  average  profiler  high 

resolution  gate  height  was  also  calculated  with  finite  difference 

derivatives  from  both  temperature  and  humidity  digitized  profiles. 

Smoothing  was  performed  as  before.  This  calculation  is  given  in  Table 

2 

7.  An  arithmetically  averaged  X  profile  from  Table  7  is  given  in 
Table  8.  High  humidity  in  the  boundary  layer  gives  extremely  high 


Table  6. 

vertical  profiles  for  each  flight  during  which  the  profiler 
operated,  calculated  at  each  aircraft  observation  level  in  MSL 
altitude  with  aircraft -measured  temperature  and  specific  humidity 
profiles . 


Flight  2  Descent  Flight  3  Descent 
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Table  6. 
(continued) 

Flight  6 
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x- 
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x- 
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Table  7. 
(continued) 
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2 

values  of  C  -radar  at  these  levels.  Aircraft  profiles  of 
n  r 

2  2  2 
-optical  adjusted  with  X  to  give  -radar  show 

2 

agreement  with  the  profiler  derived  -radar  profiles  (Figure 
26). 

For  the  two  night  flights  during  which  multiple  instruments  took 
2 

data,  comparisons  of  -optical  vertical  profiles  are  shown.  The 
2 

level  at  which  C  was  measured  by  the  AFGL  scintillometer  varied 
n  J 

because  the  instrument  was  measuring  at  different  zenith  angles.  To 

aid  in  graphical  interpretation,  for  each  night,  all  the  AFGL 

2 

scintillometer  profile  Cn  values  were  log  averaged  and  the 

heights  arithmetically  averaged,  at  each  level.  RADC  data  (Figure  14) 

was  available  for  one  night  only  and  was  omitted  for  clarity.  Figure 

2 

27  shows  a  comparison  of  -optical  from  the  averaged  AFGL 

scintillometer  profiles,  the  aircraft,  the  VanZandt  et  al.  model 

predictions,  the  AFGL  thermosonde ,  and  the  profiler  (unadjusted  raw 
2 

Cn  -radar) .  The  ascent  profile  for  the  aircraft  on  flight  6  is 
thought  to  have  been  an  incorrect  gain  setting  which  leads  to  the 
unrealistically  high  values  in  the  lower  part  of  the  profile.  Figure 

28  is  the  same  comparison  except  that  the  radar  profile  has  been 

2  2 

converted  to  -optical  using  the  calculated  X  profiles  of 

Table  7.  These  plots  show  again  the  agreement  between  the  model  and 

the  aircraft,  while  the  thermosonde  has  good  agreement  with  the 

2 

scintillometer.  The  profiler  -optical  values  (from  the  hourly 
average  observation  during  the  aircraft  flight)  show  some  agreement 


with  the  thermosonde. 


ALTITUDE  (km) 


Figure  27.  (continued) 


7.4  Propagation  Parameters 


Other  optical  parameters  of  interest  are  the  coherence  length 

rQ  and  isoplanatic  angle  6  .  These  are  related  to  vertical 
o 

integrals  of  Cn  (Murphy  and  Battles  1986) 


kr  [  C2(z)  dz  f3/5 

o 


:(z)  z5/3dz  ]'3/5 


(7.3) 


(7.4) 


where  kt  and  kf  are  functions  of  wavelength  (a  ‘  dependence). 

Assuming  a  wavelength,  and  a  linear  relation  between  observations  of 
2 

Cn  -optical  at  each  level,  the  integrals  can  be  evaluated 

numerically.  But  contributions  from  above  and  below  the  measured 

profile  must  be  calculated  based  on  theoretical  assumptions.  In  the 

2 

boundary  layer,  Murphy  and  Battles  (1986)  assumed  Cn  to  go  as 


Cn<z>  -  Cn<zb>  <Vzb>'2/3 


(7.5) 


where  z^  is  the  lowest  level  of  the  measured  profile.  In  the  upper 

2 

atmosphere,  was  modeled  by  an  exponential  decrease 


C2n(z)  -  cj(zt>  e-a(z-zt> 


(7.6) 


where  a  -  (In  10)  DR,  and  DR  (the  atmospheric  drop  off  rate)  -  1.3  X 
-4 

10  .  With  these  relationships,  the  integral  portions  for 

above  and  below  become  (Murphy  and  Battles  1986) 


Above  -  z3/3  (C3(zt)/a)  [1  +  (5/3)  (az^'1  +  (10/9)  (a2  z2  j 

(7.7) 


.V."AW?V 


Below  - 


2  Zb  Cn<zb> 


(7.8) 


and  for  r 


Above  -  Cz(zc)/a 


Below  -  3  zb  C^(Zb)  [1-  (6/zb)1/3; 


(7.9) 


(7.10) 


As  part  of  the  AFGL  thermosonde  flights,  rQ  and  8q  were 

calculated  using  equations  7.3  and  7.4  linearly  integrated  for  just 

the  thermosonde  profile,  assuming  A  -  . 3  p  (Robert  Beland,  letter  to 

William  Syrett,  September  1986;  Robert  Beland,  telephone  conversation. 

November  1987).  These  values  are  shown  in  Table  9.  r  and  8 

oo 

were  also  calculated  by  Murphy  and  Battles  (1986)  with  all  the  above 
integral  equations  and  theoretical  assumptions  for  the  atmosphere 
above  and  below,  for  the  AFGL  scintillometer  data,  using  A=-  0.5  \x .  An 
average  value  for  each  night  is  given  in  Table  10.  No  rQ  values  are 
shown  for  the  AFGL  scintillometer  because  the  boundary  layer  has  a 
large  contribution  to  the  rQ  integral.  The  lowest  scintillometer 
measurement  level  is  roughly  2  km  AGL.  Comparisons  of  contributions 
to  the  total  integral  showed  that  the  theoretical  extrapolations  below 
the  measured  profile  were  contributing  the  major  percentage  to  the 
computed  value.  Calculated  rQ  values  were  thus  considered 
unreliable  (Murphy  and  Battles  1986).  To  facilitate  comparison,  this 
same  method  and  assumptions  were  used  to  compute  rQ  and  9  from 
the  aircraft  data  at  0.3  p  and  0.5  n  wavelength  (Table  11).  Both  the 
aircraft  and  thermosonde  profiles  start  at  much  lower  AGL  altitudes 
and  thus  the  measured  profile  is  a  significant  contributor  to  the  r 
calculations  for  both  these  instruments  (Robert  Beland,  telephone 


Table  10. 

Average  of  AFGL-calculated  optical  parameters  from  scintillometer 
data . 


AFGL  Scintillometer  Data 
(calculated  with  A-=.5 n) 


GMT  Time 


5/4/86  0342-0513 


SQ  ad) 


5/6/86  0317-0420  10.29 


a=  Scintillometer  Measurement 
b-  Atmosphere  Above 
c-  Atmosphere  Below 


Percent  Contribution 


I 


Table  11. 


Optical  parameters  from  aircraft  data  calculated  with  the  Murphv  air 
Battles  assumptions. 


Aircraft  Data 
(calculated  with  A=-.5ji) 


Date 

Flieht 

0O  (^rad) 

Percent  Contribution 

4  May 

6 

8.3 

a  84.8 
b  15 
c  0.2 

6  May 

8 

7.9 

Aircraft  Data 
(calculated  with  X— 

a  49.2 
b  50.7 
c  0 

.3/0 

Date 

Flieht 

9q  (Miad) 

Percent  Contribution 

4  May 

6 

4.5 

a  84.8 
b  15 
c  0.2 

6  May 

8 

4.3 

a  49.2 
b  50.7 
c  0 

Date 

Flieht 

(cm) 

Percent  Contribution 

4  May 

6 

2.4 

a  68.5 
b  0.5 
c  31 

6  Mav 

8 

4.8 

a  82.2 

b  6 . 4 
c  11.3 


a=  Thermosonde  Measurement 
b=  Atmosphere  Above 
c=  Atmosphere  Below 
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conversation,  November  1987).  The  8q  integral,  inversely,  weights 

the  upper  atmosphere  most  heavily.  Since  the  aircraft  is  limited  in 

altitude  to  =-10  km,  the  contribution  to  8  is  limited.  However,  due 

o 

2 

to  the  integral  approximation,  profiles  in  which  Cn  drops  off 
sharply  with  altitude  will  weight  the  measured  layer  more  heavily. 

This  is  seen  in  the  profile  of  flight  6.  6  values  from  the 

scintillometer  and  the  aircraft  (at  A  -  0.5  n)  were  similar  on  4  May, 
but  disagreed  on  6  May.  The  thermosonde  and  aircraft  8  values  (at 
A  -  0.3  /j)  disagreed  on  4  May,  but  agreed  on  6  May. 

rQ  was  also  measured  at  optical  wavelengths  during  the  EWAK 
experiment  by  the  NPS  rQ  scintillometer.  The  measured  rQ  average 
for  the  aircraft  flight  time  on  the  night  of  4  May  1986  was  5  cm.  The 
thermosonde  and  aircraft  rQ  values  (at  A  -  0.5  /i)  differed  by  a 
factor  of  two  on  4  May.  The  night  of  6  May,  the  thermosonde  and 
aircraft  rQ  values  were  similar. 

7.5  VanZandt  et  al .  Model  for  e 

The  VanZandt  et  al .  model  €  predictions  are  compared  to  e 

profiles  from  the  aircraft  data  (Figure  29).  The  hot  wire  failed 

during  flight  4,  so  that  only  a  small  amount  of  data  is  available. 

The  model  shows  predictions  generally  ranging  in  the  area  of  5  M  10-- 

for  all  flights.  There  is  poor  agreement  between  the  aircraft  and 

model,  the  aircraft  (except  for  flight  3)  often  giving  values  three  or 

four  orders  of  magnitude  smaller.  Aircraft  values  also  varv  much 

-  3  -  7 

more,  ranging  from  10  to  10 

A  crude  average  profile  of  r\ ,  the  microscale,  based  on  the 

3  25 

aircraft  e  profiles  is  shown  in  Figure  30.  r)  -  {v  /e)' 


was 


computed  for  each  e  value  using  the  previously  mentioned  formula  for 


with  the  standard  atmosphere  temperature  profile,  and  then  averaged 
1  km  layers.  The  average  value  of  r?  spans  one  order  of  magnitude, 
from  1  mm  in  the  boundary  layer  to  an  average  of  1  cm  above  the 
boundary  layer . 


CHAPTER  8 


SUMMARY  AND  CONCLUSIONS 


The  atmosphere  has  a  temporally  and  spatially  variable  refractive 

index  due  to  its  inherently  turbulent  nature.  This  can  affect 

electromagnetic  beam  propagation  by  introducing  beam  steering,  image 

dancing,  beam  spread,  spatial  coherence  degradation,  temporal 

coherence  degradation,  and  scintillation  (Dewan  1980) .  Refractivity 

is  wavelength  dependent.  At  microwave  radar  wavelengths,  molecular 

composition  (e.g.,  water  content)  dominates  refractivity  changes.  At 

visible  wavelengths,  density  (parameterized  by  temperature)  dominates 

(Balsley  and  Gage  1980).  For  many  quantification  purposes,  the 

2 

refractive  index  turbulent  structure  constant,  C  ,  is  a  key 

parameter.  This  parameter  can  be  inferred  from  optical  turbulence 

(scintillometer) ,  temperature  turbulence  (aircraft  instrumentation  or 

thermosonde)  or  radar  backscatter  (profiler)  measurements.  However, 

2 

Cn  from  a  radar  profiler  will  differ  because  of  the  wavelength 

dependence  of  refractivity.  Theoretically,  this  can  be  compensated 

for,  if  the  humidity  profile  is  known  (VanZandt  et  al.  1981). 

2 

Vertical  integrals  of  give  transverse  coherence  length,  r  . 
and  isoplanatic  angle,  SQ.  These  can  also  be  measured  remotely 
(based  on  optical  turbulence)  by  an  rQ  scintillometer  and 
isoplanometer ,  respectively. 

During  an  atmospheric  optics/meteorology  experiment  (acronym 
EWAK)  conducted  at  Penn  State  University  primarily  during  April  and 
May  of  1986,  data  was  collected  by  the  aforementioned  instruments. 
AFGL  and  RADC  scintillometers  produced  vertical  profiles  of 


2 

-optical.  An  AFGL  thermosonde  produced  vertical  profiles  of 

2 

C^,  and  other  meteorological  variables.  A  PSU  profiler  produced 

2 

vertical  profiles  of  wind  direction  and  speed,  and  -radar.  The 

instrumented  ARA  research  aircraft  produced  vertical  profiles  of 

2  2 

temperature  and  velocity  turbulence  (C^,  and  )  and  other 
meteorological  variables. 

An  indepth  analysis  of  the  aircraft  turbulence  data  was 

performed.  The  turbulence  instrumentation  aboard  consisted  of 

cold-wire  and  hot-wire  sensors  (for  temperature  and  velocity 

variations ,  respectively)  and  FM  recording  apparatus .  The  taped  data 

was  processed  via  FFT  to  produce  one -dimensional  variance  spectra 

(wavenumber  range  0.01  to  10  m'^).  Flights  usually  produced  a  10  km 

vertical  profile;  the  data  was  processed  to  give  roughly  0 . 5  km 

vertical  resolution,  similar  to  that  of  the  profiler. 

An  atmospheric  filter  function  was  derived,  based  on  PBL  spectra 

and  -5/3  inertial  subrange  slope.  This  function  was  applied  to 

compensate  for  low  frequency  dropoff,  as  part  of  a  spectral  editing 

program.  Editing  was  based  on  percent  error  between  a  regression 

analysis  and  theoretical  -5/3  slope.  Noise  subtraction  was  performed 

in  some  low  signal  cases.  In  general,  PBL  spectra  clearly  showed  the 

classical  -5/3  inertial  subrange  slope,  out  to  the  broadband 

instrument  noise  level.  A  majority  of  spectra  in  the  free  tropospher 

showed  good  evidence  of  the  inertial  subrange  with  -5/3  slope. 

Further  data  collection  should  emphasize  the  low  frequencies,  since 

this  appears  to  be  the  low  noise  region.  Occasional  anomalous  and  lo 

signal-to-noise  spectra  were  encountered  above  the  PBL.  Low  signal 

2  2 

spectra  represent  quiescent  layers .  C_  and  C  were 


calculated  from  the  regression  fit  to  the  inertial  subrange  power 
spectral  density.  The  rate  of  dissipation  of  turbulent  kinetic 

energy,  e,  can  be  calculated  by  using  the  so-called  Corrsin  relation 

2  2 
for  velocity  turbulence.  C  can  be  calculated  from  C„  . 

J  n  T 

Considerable  interest  has  developed  in  models  that  relate 

microturbulence  parameters  to  the  mean  gradients.  One  such  model, 

proposed  by  VanZandt  et  al .  (1978;  1981),  has  as  a  key  variable  the 

gradient  Richardson  number,  R^.  In  actively  turbulent  regions, 

2  2 

can  be  related  to  C^,  and  as  follows: 

4  /  cl  -  1.6  Rj/  <Pr  -  R.)  (8/ g)  (38/3z)  (8.1) 

2  2 

A  scatter  plot  of  CT  versus  (8/ g  38/ dz)  in  the 

free  troposphere  has  a  slope  that  corresponds  to  1.6  R./(P  - 

R^).  This  plot  showed  a  range  of  values  from  0.3  to  10  for  high  and 

2  2 
low  values  of  CT  ,  respectively.  Scatter  plots  of  CT  versus 

2 

Cu  showed  high  correlation  of  these  two  parameters ,  with  log 

2  2 

Cu  -  1.7  log  C.j,  .  This  implies  that  the  mean  gradient 

structure  of  the  free  troposphere  would  have  a  probability  densitv 

maximum  (Fairall  and  Markson  1985).  This  could  be  illustrated  as  a 

single  curve  in  N-S  space,  where  N  is  the  Brunt- Vaisala  frequency 

squared,  and  S  is  the  shear  squared  and  R^  -  N/S .  Plots  of  this 

type  were  done  using  aircraft  temperature  data  and  profiler  wind  dat 

A  distribution  was  obtained  roughly  corresponding  to  the  curve 

depicted  by  Fairall  and  Markson  (1985).  Generally,  regions  of  high 

shear  are  at  or  below  the  0.25  line,  as  expected. 

2  2 

Plots  of  C.j,  /  (8/ g  38/3 z)  versus  a  turbulent  activity 

parameter  (Gregg  1987)  clearly  show  agreement  with  the  suggested 


activity  levels  and  associated  values.  As  the  ordinate  value 
[corresponding  to  1.6  R^/(P  -  R^)]  decreases,  turbulent 

activity  increases.  At  high  activity  levels,  the  ordinate  value 
approaches  0.4  (equivalent  to  1.6  R^  when  approximating  Pr  as  1 
and  R.^  as  0.25). 

2 

Values  of  N  and  S  were  used  to  obtain  values  of  C  and  e 

n 

2 

with  the  model  of  VanZandt  et  al .  Comparisons  of  Cn  profiles 
measured  by  the  various  instruments  showed  good  agreement  between 
scintillometer  and  thermosonde,  and  between  the  VanZandt  et  al.  model 

and  the  aircraft  data.  The  profiler  (converted  optical  values)  had 

2 

some  agreement  with  the  thermosonde.  The  difference  in  Cn 

profiles  between  methods  varied.  An  average  ratio  of  each  profile 

versus  the  thermosonde  profile  gave  the  following  factors  for  flight 

6:  aircraft  11.7  (standard  deviation  23.4),  model  11.5  (stand,  dev. 

18.5),  profiler  4.3  (stand,  dev.  5.3),  and  the  scintillometer  1.6 

(stand,  dev.  0.5).  In  the  lower  levels  of  the  flight  6  ascent,  the 

aircraft  profile  is  thought  to  have  an  error  in  gain  setting.  For 

flight  8,  the  factors  were:  aircraft  4.2  (stand,  dev.  4.5),  model  2. 

(stand,  dev.  2),  profiler  12  (stand,  dev.  34),  and  the  scintillometer 

0.7  (stand,  dev.  0.6).  Similar  differences  between  aircraft  and 

thermosonde  were  noted  by  Brown  and  Good  (1984).  It  was  later  though 

that  this  difference  might  have  been  due  to  spectral  analysis 

technique.  However,  the  increased  effort  to  improve  spectral  analvs: 

in  this  paper  apparently  seemed  to  have  little  effect.  Structure 

constant  values  obtained  by  the  two  spectral  analysis  methods  examine 

in  this  work  (FFT  board  and  the  signal  analyzer)  differed  on  the 

2 

average  by  about  a  factor  of  two  for  C_  ;  however  individual 


A 


values  differed  by  as  much  as  two  orders  of  magnitude . 

2 

To  account  for  the  difference  in  profiles  between 

instrument  systems,  further  research  might  seek  to  determine  the 

difference  in  values  obtained  by  thermosonde  and  aircraft 

instrumentation,  since  they  actually  measure  some  of  the  same 
2 

variables  (e.g.,  )  in  slightly  different  ways.  Since  the 

scintillometer  tended  to  agree  with  the  thermosonde  and  the  VanZandt 
et  al .  model  tended  to  agree  with  the  aircraft,  no  clear  choice  is 
evident.  The  radar  does  not  have  an  absolute  calibration,  so  there  is 
no  significance  to  its  agreement  with  either  system  (Christopher 
Fairall,  personal  communication,  December  1987). 

Another  noticeable  difference  between  the  aircraft  and 
thermosonde  is  the  measured  vertical  profiles  of  potential 
temperature.  Since  potential  temperature  and  its  gradient  are  input 
parameters  in  the  VanZandt  et  al .  model,  this  difference  could 
obviously  affect  model  output.  Barat  and  Bertin  (1984)  state  the  need 
for  accurate  temperature  profiles  and  suggest  changes  in  Richardson 
number  depend  more  on  temperature  gradient  than  on  shear.  This  study 
only  used  aircraft  data  as  model  input.  Comparing  model  output  from 
thermosonde  and  aircraft  input  could  be  another  interesting  comparison 
for  further  study.  There  is  also  an  obvious  need  for  additional 
simultaneous  data  collection.  Logistical  and  weather  problems 
contributed  to  make  only  two  nights  of  simultaneous  operation  possible- 
during  EVAK. 

2 

Numerical  integration  of  C  gave  estimates  of  r  and  9 
&  n  °  o  o 

Values  from  scintillometer,  thermosonde  and  aircraft  did  not  show 
consistent  agreement.  Based  on  the  aircraft  values,  both  the 
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scintillometer  and  thermosonde  differed  by  an  average  of  17%  for  6  , 

and  the  thermosonde  differed  by  an  average  of  35%  for  r  .  With  such 

a  limited  data  set  it  is  not  possible  at  this  time  to  conclude  whether 

2 

profiles  of  can  be  used  to  estimate  these  optical  parameters. 
Acceptable  approximations  for  unmeasured  segments  of  the  atmosphere 
must  also  be  further  studied. 

2 

Despite  good  agreement  on  ,  the  aircraft  data  and  the 
VanZandt  et  al .  model  output  clearly  disagreed  on  vertical  profiles  of 
e.  This  difference  was  also  noted  by  Fairall  and  Markson  (1985). 


They  suggest  this  effect  is  due  to  the  adjustment  of  the  model 

parameter  L,  turbulent  layer  thickness,  to  produce  the  best  fit  of 

2 

model  output  to  profiler  data.  Since  few  measurements  of  e 

have  been  done,  there  is  a  definite  need  for  a  larger  data  base. 

Other  areas  of  study  might  include  adjusting  L  to  best  fit  aircraft  « 

profiles.  However,  this  would  require  either  adjusting  the  constants 

(a  and  b)  in  the  model  equations  or  modifying  the  theory  in  order  to 

2 

preserve  the  good  predictions  for  .  Raw  high  speed  aircraft 
data  could  also  be  used  to  study  turbulent  "episodes"  individually 
without  vertical  averaging.  The  fraction  of  the  profile  that  is 
actually  turbulent  could  be  computed  and  the  microturbulence  ratio 
could  be  examined  for  each  active  layer. 

Overall,  EWAK  provided  a  chance  to  compare  methods  of  measuring 

2 

,  as  well  as  an  opportunity  to  obtain  microturbulence  data  that 
has  some  relation  to  the  mean  gradient  structure.  Hopefully,  further 
research  will  provide  a  relationship  that  could  be  exploited  in  fur-:.- 


predictive  models. 
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